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ABSTRACT 


Double-exposure  specklegrams  were  made  of  aluminum 
2024-T4  rectangular  plates.  The  double -exposure  speckle 
technique  was  applied  to  the  aluminum  plates  to  measure 
in-plane  translation  and  rotation  and  out-of-plane  tilt. 

Strain  and  vibrational  analyses  were  also  conducted.  The 
vibrational  results  were  compared  with  results  previously 
obtained  using  double -exposure  time  average  holography  to 
ascertain  their  applicability  to  Naval  Engineering. 

Background  in  Fourier  transforms,  wave  theory,  Fraunhofer 
diffraction  theory,  and  other  areas  of  physical  optics  were 
investigated  to  enhance  the  development  of  speckle  metrology 
techniques. 

Results  of  the  experiments  indicate  that  speckle  metrology 
used  in  conjunction  with  holography  has  significant  ootential 
for  Naval  Engineering  apolications.  Further  research  is  re¬ 
quired  to  develop  these  techniques.  These  procedures  currently 
are  being  utilized  extensively  in  the  aerospace  industry  and 
automotive  research  and  development  areas. 
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DEFINITION 


a  Angular  measurement 
diverging  core 


NOMENCLATURE 


UNITS 


degrees 


and  normal  to 
(See  Fig.  4) 

Phase  difference  be* 
tween  2  waves 
Phase  difference 
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equation  numbers 
139 


a 

Angle  of  prism  wedge 

degrees 

141,  143 

Y 

,  Angle  of  tilt 

degrees 

128,  136,  137,  138 

Yx 

y Shear  strain 

- 

145 

6 

Delta  function 

- 

47,  49,  50,  51,  52, 

53,  56,  57,  58,  59, 

60,  61,  62 

6 

Angle  between  P  P 
o 

degrees 

19,  24 

100,  101,  107 


144,  145 


r  A 

Lens  thickness 

meter 

76, 

77, 

78, 

79, 

81, 

83 

e 

Strain 

- 

145 

c 

!* 

Electrical  permittivity 

- 

5 

e 

Radius  of  integration 

meter 

12 

i'  c 

Rectangular  coordinate 

- 

20, 

22, 

23, 

24, 

25, 

27, 

28, 

29, 

30, 

32, 

' 

33, 

34 

"1 


k 


Phase  angle 

degrees 

98, 

101, 

119 

»,  120, 

129 

,  132 

n 

Axis  direction 

- 

20, 

21, 

22, 

23,  24, 

25, 

26, 

27, 

28,  29, 

30, 

32, 

33 

\ 

Wave  length 

meter 

18, 

19, 

24, 

29,  30, 

32, 

44, 

59, 

60,  61, 

87, 

88, 

89, 

90,  91, 

92, 

93, 

94, 

95,  110 

114, 

122, 

123, 

124 

125, 

126, 

128, 

129 

132, 

136, 

138, 

140 

144 

141,  143 

I 

1/second  2,  59,  60,  61 

36,  38,  40 

1/second  1,  2,  4,  72,  73,  118, 

119,  120 

«  Circular  coordinate  -  37,  38,  40,  43,  44 


V  Index  of  refraction 

of  light 

v  Temporal  frequency 

p  Axis  coordinate  for 

polar  coordinates 

(d  Circular  frequency 


ft 


axis 

Element  of  solid  angle 
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I.  INTRODUCTION 


A.  BACKGROUND 

Laser  Speckle,  that  granular,  often  annoying  pattern  pro¬ 
duced  upon  laser  illumination  of  a  diffusely  reflecting 
object,  came  into  prominence  as  a  holographer ' s  nightmare. 
While  continuing  in  that  role,  it  is  recognized  that  laser 
speckle,  used  in  double-exposure  form,  is  extremely  useful 
for  speckle  photography.  Laser  speckle  is  perhaps  the 
easiest  form  of  laser  metrology  to  exploit  because  in  its 
simplest  form,  the  equipment  required  for  speckle  image 
forming  and  processing  requires  only  single  beam, 
double-exposure,  which  eliminates  the  need  for  vibrational 
isolation.  Techniques  for  reduction  of  the  specklegram  have 
been  developed  successfully  and  applied  to  the  measurement 
of  surface  roughness  /Tl,  18,  21 7,  vibration  /S’,  20,  30,  51, 
537  and  deformation  /T2,  537,  and  displacement  /"*4,  7,  10, 

15,  38,  537  and  strain  /”3,  4,  9,  15,  31-34,  40,  537.  A  most 
convenient  form  of  data  reduction  for  speckle.  photographs  is 
based  upon  Thomas  Young's  demonstration  in  the  early  1800's 
of  the  interference  of  light  emerging  from  a  dual  pinhole 
arrangement  ^”13 7.  Thus  a  new  field  of  speckle  metrology 
evolved. 

B.  THESIS  OBJECTIVE 

The  primary  objectives  of  this  thesis  were  to  develop 
various  techniques  for  measuring  in-plane  and  out-of-plane 


displacement  and  analysis  of  strain  and  vibrations  using 
speckle  metrology.  The  objective  of  using  speckle  techniques 
instead  of  holography  procedures  in  Naval  Engineering  is  to 
enhance  the  applicability  of  photographic  metrology  to  the 
industrial  environment  /“ 16,  21,  44,  557*  Holography,  as  a 
form  of  engineering  stress,  strain,  and  vibrational  metrology, 
is  excellent  but  requires  highly  sophisticated  equipment, 
close  tolerance  set  up  requirements,  and  a  vibrationless 
laboratory  atmosphere  for  exact  measurement  /” 1,  8,  21,  44, 
537.  With  the  development  of  speckle  photography,  a  single 
beam,  easy  to  setup,  high  tolerance  system  of  engineering 
metrology  is  capable  of  being  achieved  in  an  industrial  en¬ 
vironment.  /"" 21,  4  £7 


II.  THEORY 


A.  INTRODUCTION  TO  SPECKLE  METROLOGY 

Speckle  metrology  spans  a  rather  wide  variety  of  tech¬ 
niques,  all  of  which  generally  fall  into  one  of  two  categories: 
direct  laser  (speckle)  photography  or  speckle  interferometry. 
/""2,  21,  4 47  The  distinction  between  the  two  is  somewhat 
subtle  since  both  involve  photography  and  both  involve  inter¬ 
ference.  That  is,  speckle  photography  depends  on  the  optical 
properties  of  the  speckle  which  involves  a  condition  of  inter¬ 
ference.  After  orocessing  the  speckle  photograph,  the 
analysis  of  the  interference  yields  fringes.  These  fringes 
are  directly  correlated  to  the  object  motion.  Therefore, 
there  is  no  great  difference  between  speckle  photography  and 
speckle  interferometry.  /— 21,  4 £7  Hence,  if  there  are  regions 
in  the  two  images  where  speckle  patterns  of  each  are  well 

J 

correlated,  let  the  process  be  referred  to  as  speckle 
photography.  If,  instead,  the  fringes  form  as  a  result  of 
fluctuations  in  the  correlation  of  the  speckle  patterns  be¬ 
tween  two  images,  whether  or  not  there  is  translation  be¬ 
tween  the  correlated  portions  of  the  pattern,  let  the  process 
be  called  speckle  interferometry. 

B.  ELECTROMAGNETIC  WAVE  THEORY 

Light,  being  a  form  of  electromagnetic  radiation,  is 
characterized  by  its  amplitude,  wave  length  (or  frequency), 
phase,  polarization,  speed  of  propagation,  and  direction  of 


26 
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propagation.  /—13,  22,  43,  537  Light  can  be  described  by 
specifying  the  temporal  and  spatial  dependence  of  its  electric 
field  vector,  E.  The  simplest  type  of  electromagnetic  wave 
is  the  linearly  polarized  plane  wave  which  has  been  assumed.  . 
If  such  a  wave  is  polarized  in  the  y  direction  and  propagates 
in  the  z  direction,  the  three  components  of  E  are: 

Ez  -  0, 

Ey  *  A  cos  (ut  -  kz),  l  * ) 

E,  -0. 


where  A  is  the  amplitude  of  the  wave. 

w  and  wave  number  K  are  given  by 

a  »  2irr 


k 


2sr 
X  ’ 


The  circular  frequency 


(2) 

(3) 


where  v  is  the  temnoral  frequency  and  X  is  the  wave 
length.  The  frequency  of  light  is  on  the  order  of  1015HZ. 

More  generally,  the  direction  light  travels  is  described 
by  its  propagation  vector  K  which  has  a  magnitude  K  =  2ir/X 
and  points  in  the  direction  of  propagation.  A  plane  wave  is 
a  wave  whose  phase  at  any  instant  of  time  is  constant  at  all 
points  normal  to  K.  If  r  ■  xi  +  yT  +  xic  is  the  position 
vector  of  any  point  in  space,  as  shown  in  Fig.  1,  the 
equation  of  a  linearly  polarized  plane  wave  is: 


-  0. 

Ey'  -  *4cos(uf-kf)> 

£,.-0. 


(4) 


A  surface  over  which  phase  is  constant,  in  this  case  the 


planes  K»r  =  CONSTANT,  is  called  a  wave  front. 


Figure  1.  Plane  wave  of  light  propagating  in  the  direction 
specified  by  the  propagation  vector  K.  From  Vest 

/~537 

Lasers  emit  narrow  beams  of  monochromatic  light  with 
almost  perfectly  plane  wave  fronts.  Most  lasers  emit 
linearly  polarized  light.  Light  from  a  typical  He-Ne  con¬ 
tinuous  wave  (CW)  laser  has  a  wavelength,  X  *  632.8  nm.  The 
laser  emits  a  beam  of  about  5  mm.  in  diameter.  Laser  liaht, 
as  was  mentioned  previously,  is  highly  coherent. 

It  was  also  noted  previously  that  the  frequency  of  light 
was  on  the  order  of  1015Hz.  Practical  detectors,  such  as  the 
eye  and  photographic  film,  are  not  capable  of  responding  to 
such  extremely  rapid  variations.  Rather,  they  respond  to 
irradiance  I  which  is  the  time-average  energy  flux  of  the 
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light  wave.  Using  electromagnetic  wave  theory,  it  can  be 
shown  that 


/  -  eu<E2>, 


(5) 


where  e  is  the  electrical  permittivity  of  the  medium  in 
which  light  travels;  v  is  the  speed  of  propagation.  The 
key  point  is  that  the  irradiance  is  proportional  to  the  tine 
average  of  the  electrical  intensity  vector  squared;  the  time 
average  is  indicated  by  <>  . 

C.  DIFFRACTION  THEORY 

1.  Derivation  of  the  Fresnel-Kirchhof f  Diffraction 
Formula' 

According  to  the  HUYGEN-FRESNEL  theory  /T 37,  the  light 
disturbance  at  a  point  P  arises  from  the  superposition  of 
secondary  waves  that  proceed  from  a  surface  situated  between 
this  point  and  the  light  source. 

Consider  a  strictly  monochromatic  wave 

V(x,y,z,t)  «  U(x,y,z)  e”iwt  (6) 

In  a  vacuum,  the  space-dependent  part  satisfies  the 
time-independent  (Homogeneous)  wave  equation. 

(V2  +  k2)  U  »  0  (7) 

where  K  ■«/ c.  This  is  known  as  the  Helmholtz  equation. 

Let  v  be  the  volume  bounded  by  a  closed  surface  S, 
and  let  P  be  any  point  within  the  surface;  (see  Fig.  2), 
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Assume  that  U  possesses  continuous  first  and  second  order 
partial  derivatives  within  and  on  the  surface.  If  U'  is 
any  other  function  which  satisfies  the  same  continuity 
requirements  as  U,  then  by  Green's  theorem 

JJJVw -mu,.  -JX(Fiir-  v'%)  m 


where  3()/3n  denotes  differentiation  along  the  inward  normal 
to  S.  If  U*  also  is  a  solution  to  the  Helmholtz  equation, 
i.e.,  if 


(V*  +  V)ir  =»  o, 


(9) 


then  it  follows  that  the  integrand  on  the  left  of  equation 
(8)  vanishes  at  every  point  of  v  and  consequently 


dS  =  0. 


(10) 


Figure  2.  Derivation  of  the  Helmholtz-Kirchhoff  integral 
theorem:  region  of  integration.  From  Born  and 
Wolf  /“13  7 
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Let  U* (x,y,z)  =  e  /s,  which  is  a  spherical  wave 
converging  to  the  point  P(s*0),  where  s  denotes  the  distance 
from  P  to  the  point  (x,y,z).  This  function  has  a  singular¬ 
ity  at  s  *  0;  and  since  U1  was  assumed  to  be  continuous  and 
differentiable#  P  must  be  excluded  from  the  domain  of  in¬ 
tegration.  Surrounding  P  with  a  small  sphere  of  radius  e,  the 
integration  is  expected  throughout  the  volume  between  S  and 
the  surface  S'  of  this  sphere.  In  place  of  equation  (10) 


(11) 


whence 


,12) 


where  A  denotes  an  element  of  the  solid  angle.  Since  the 

i 

integral  over  S  is  independent  of  e ,  the  integral  on  the 
right-hand  side  may  be  replaced  by  its  limiting  value  of 
c  ♦  o.  The  first  and  third  terms  in  the  integral  give  no 
contribution  in  the  limit,  and  the  total  contribution  of 
the  second  term  is 


V(P) 


(13) 


This  equation  is  known  as  the  Kirchhoff  integral 
theorem.  It  relates  the  value  of  any  scalar  wave  function 
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at  any  point  P  inside  an  arbitrary  closed  surface  to  the 


value  of  the  wave  function  at  the  surface. 

Note  that  as  K-*0,  the  time-independent  wave 

2 

equation  (9)  reduces  to  Laplace's  equation  7  U  *  0,  and 
equation  (13)  goes  into  the  formula  of  potential  theory. 


U(P) 


(14) 


Only  monochromatic  waves  will  be  analyzed,  but  a 
general  form  of  Kirchhoff's  theorem  could  be  derived  which 
applies  to  waves  not  necessarily  monochromatic. 

Consider  a  monochromatic  wave  from  a  point  source 
PQ  propagated  through  an  opening  in  a  plane  opaque  screen, 
and  let  P  be  the  point  at  which  the  light  disturbance  is 
to  be  determined.  It  is  assumed  that  the  linear  dimensions 
of  the  opening,  while  large  compared  to  the  wave  length  of 
light,  are  small  compared  to  the  distances  of  both  PQ  and  P 
from  the  screen. 

To  find  the  disturbance  at  P,  Kirchhoff's  integral 
is  taken  over  the  surface  S  formed,  as  seen  in  Fig.  3,  by: 
(1)  the  opening  A,  (2)  a  portion  B  of  the  non-illuminated 
side  of  the  screen,  and  (3)  a  portion  C  of  a  large  sphere  of 
radius  R,  centered  at  P  which,  together  with  A  and  B,  forms 
a  closed  surface. 

Kirchhoff's  theorem  then  gives 


where  S  is  the  distance  of  the  element  ds  from  P  and 


denotes  the  differentiation  along  the  inward  normal  to 
the  surface  of  integration. 

The  values  of  U  and  3U/3n  on  A,  B,  and  C  which 
should  be  substituted  into  equation  (15)  are  never  exactly 
known.  However,  it  is  assumed  that  everywhere  on  A,  except 
in  the  neighborhood  of  the  rim  of  the  opening,  U  and  3U/3n 
will  not  appreciably  differ  from  the  values  obtained  in  the 
absence  of  the  screen  and  that  on  B  these  quantities  will 
be  equal  to  zero. 


Then  Kirchhoff  set 

on  jS:  U  =  U">, 

on  38  \  V  =  0 


3£  dUw 
dn  dn 


(16) 


where 

<17 

are  the  values  of  U  relating  to  the  incident  field  (see 
Fig.  3b)  and  A  is  a  constant.  The  approximations  of  equa¬ 
tion  (16)  are  called  Kirchhoff1 s  boundary  conditions  and 


are  the  basis  of  Kirchhoff's  diffraction  theory. 


<a> 


(a)  Surface  of  integration  (b)  Definition  of  angles 

Figure  3.  Illustrating  the  derivation  of  the  Fresnel- 

Kirchhoff  diffraction  formula.  From  Born  and 
Wolf  /“13 J 

It  remains  to  consider  the  contribution  from  the 
spherical  portion  C.  If  the  radius  R  is  sufficiently  large 
so  that  at  the  time  when  the  disturbance  at  P  is  considered 
no  contributions  from  C  could  have  reached  P,  then  the 
integral  over  C  will  vanish.  Neglecting  the  normal  terms 
1/r  and  1/s  in  comparison  to  k,  the  Fresnel-Kirchhoff 
diffraction  formula  is  obtained 

_  \A  fr  e“(rfrt 

u(p)  -  -  j}  J  J^TT*  t008  (n- r)  ~  008  (n*  ( 18 ) 

Closer  consideration  of  the  Fresnel-Kirchhoff  dif¬ 
fraction  integral,  equation  (18)  ,  shows,  as  the 
element  ds  explores  the  domain  of  integration,  (r  +  s)  will 
in  general  change  by  many  wave  lengths,  so  that  the  factor 


will  oscillate  rapidly.  On  the  other  hand. 


ik  (r  +  s) 
e 

the  distances  of  the  points  PQ  and  P  from  the  screen  are 
large  compared  to  the  linear  dimensions  of  the  aoerture, 
and  therefore  the  factor  /“cos  (n,r)  -  cos  (n,s)7  will 
not  vary  appreciably  over  the  aperture.  All  angles  are 
assumed  to  be  small.  The  factor  /“cos  (n,r)  -  cos  (n,s^_7 
may  be  replaced  by  cos  6  ,  where  5  is  the  angle  between  the 
line  PQP  and  the  normal  to  the  screen  (as  seen  in  Fig.  4). 
Finally,  the  term  1/rs  may  be  replaced  by  1/r's',  where 
r'  and  s'  are  the  distances  of  Pq  and  P  from  the  origin  re¬ 
spectively.  The  Kirchhoff- Fresnel  equation  then  reduces 
to 


U(P) 


-4i  cos  6  f  f 

tttJJ, 


(19) 


Cartesian  coordinates  are  used  with  the  origin  in 
the  aperture  at  0  and  with  the  x  and  y  axes  in  the  plane 
of  the  aperture  and  the  positive  z  direction  to  point  into 
the  half  space  that  contains  the  point  P  of  observation  (see 
Fig.  4). 


Figure  4.  Diffraction  at  an  aperture  in  a  plane  screen. 
Modified  from  Born  and  Wolf  £“13_7 


1 


If  (x  y  z  )  and  x,y,z)  are  the  coordinates  of  P„ 
o  o  o  o 

and  P  respectively,  and  (  5,n  )  are  the  coordinates  of  a 
point  Q  in  the  aperture,  then 


**  =  (*o  -  £)*  4-  (y«  -  »?)*  -f  V. 
**  =  (*-  #)*  4-  (y  -  »?)*  4-  **, 


(20) 


r'*  =  *,*  +  !/o*  +  *»’•  (21) 

a'*  =  z*  4-  y*  4-  2‘- 


r*  *.  r'»  —  2(z^  4-  i'o>?)  4-  £*  4-  *?*, 
**»*'*-  2(x*  4-  M)  +  P  4-  n*. 


(22) 


Since  the  linear  dimensions  of  the  aperture  are  small  con- 
pared  to  both  r'  and  s',  r  and  s  may  be  expanded  as  a  power 
series  in  £/r',  n/r' ,  5/r*  ,  n/r'. 


r  ~  r* 


4-  y«»7  ,  I*  4-  »?*  _  (*>l  4-  yp^)* 
r*  +  2r'  2r'» 

*£  4-  y>?  ,  f*  4-  >7*  (*£  4-  yn)*  _ 

a'  +  2*'* 


(23) 


Substituting  from  equation  (23)  into  equation  (19)  gives 


17(P) 


*  008  j*  /{£({>;, 


(24) 
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where 


/(*.<?)-  - 


V  +  y*n 


*£  ±  M  ,  £*  +  V*  .  £*  +  *}* 
a'  *  2r'  +  2*' 


(*o£  +  y«q)*  _  (g£  +  v*;)* 
2r'*  <  &»'* 


(25) 


The  first  two  directional  cosines  are  (1  ,  m  )  and  (l,m) 

o  o 


(26) 


and  equation  (27)  may  be  written  in  the  form 


/(£.  *l)  =»  (l,  —  l)i  +  (m0  —  m)t]  +  i  ({*  +  rj*)  _  +  "V?)*  __  (j£  +  ”»?)*  (27) 

The  problem  of  determining  the  light  disturbance  at 
P  has  been  reduced  to  the  evaluation  of  the  integral  in 
equation  (24) .  The  evaluation  can  be  simplified  for  two  cases 
of  particular  interest.  The  quadratic  and  higher  terms  in  £ 
and  n  may  be  neglected  in  f  for  Fraunhofer  diffraction?  when 
the  quadratic  terms  cannot  be  neglected,  Fresnel  diffraction 
is  being  represented.  Fortunately,  the  simpler  case  of 
Fraunhofer  diffraction  is  of  much  greater  importance  in 
optics.  /“13  7 
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2.  Fraunhofer  Diffraction 


Fraunhofer  diffraction  occurs  when  both  the  incident 
and  diffracted  waves  are  effectively  plane.  This  will  be  the 
case  when  the  distances  from  the  source  to  the  diffracting 
aperture  and  from  the  aperture  to  the  receiving  point  are 
both  large  enough  so  that  the  curvatures  of  the  incident  and 
diffracted  wave  can  be  neglected  (see  Fig.  5a). 

If  either  the  source  or  the  receiving  point  is  close 

enough  to  the  diffracting  aperture  so  that  the  curvature  of 

the  wave  is  significant,  then  Fresnel  diffraction  occurs  (see 
Fig.  5b) . 


(a)  Fraunhofer  case  (b)  Fresnel  case 

Figure  5.  Diffraction  by  an  aperture.  From  Cowley  /J7_7 
The  wave  is  effectively  plane  over  the  aperture  if 
the  second  order  and  higher  terms  are  negligible.  Strictly 
speaking,  these  terms  disappear  only  in  the  limiting  case 
r'-**,  s 1  -*■  «  ,  i.e.,  when  both  the  source  and  the  point 
of  observation  are  at  infinity.  It  is,  however,  evident 
that  the  second  order  terms  do  not  contribute  appreciably 
to  the  integral  if 
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(V  4-  m0 r})*  (l£  +  tfh;)*  I 


2jt. 


The  conditions  under  which  equation  (28)  will  be 


satisfied  exist  when 


’  \f\ >  (i*  ±-p2X  ^  |,»|  ^  (JL+p g 


or  if 


i  i  K|A 

_  +  -  =  0  and  V, m0*, l\ m*<  ^  +  ij*)^ 


The  conditions  in  equation  (29)  give  estimates  for 
r*  and  s'  for  which  the  Fraunhofer  representation  may  be 


used. 


In  the  case  of  Fraunhofer  diffraction,  the  four  quan¬ 


tities  1  ,  m  ,  l,  and  m  enter  equation  (27)  only  in  the 
o  o 

combinations 


p  w  l  —  lt,  xr  m  —  mr 


The  governing  Fraunhofer  diffraction  integral  may 


then  be  written 
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V[P)  =  cjjj  jtdt], 


?E  =  IC1*1>  <32 

where  c  is  the  constant  term  in  front  of  the  integral  of 
equation  (24). 

Fraunhofer  diffraction  patterns  for  apertures  of 
various  forms  (rectangle  and  circular,  specifically)  will 
be  studied  in  detail. 

a.  Rectangular  aperture 

Consider  first  a  rectangular  aperture  of  sides 
2a  and  2b.  With  origin  0  at  the  center  of  the  rectangle 
and  with  n  and  £  axes  parallel  to  the  sides  (see  Fig.  6) . 


7 


V - 

L - Id - * 

Figure  6.  Rectangular  aperture.  From  Born  and  Wolf  /  13  7 


The  Fraunhofer  diffraction  integral  then  becomes 


V(P)  *mc[  f  e-*P*+ti>d{drj  =  cf"  f*  e-«v» drt 

J (33) 
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with  a  similar  form  for  the  integral  of  equation  (33)  taken 
with  respect  to  dn  from  -b  to  b. 

Hence  the  intensity  is  given  by 


/sin  kpay 

/sin  kqb\ 

1  '  l  ii»  J 

[  kqb  ) 

y  =  (sin 

x/x)2 

1  maximum 

of  y  ; 

Fig.  7.  It  has  a  orincipal  maximum  of  y  s  1  at  x  =  0 
and  zero  minima  at  x  *  ±ir,  +2ir,  ±3ir...  The  roots  asymptoti- 


Figure  7.  Fraunhofer  diffraction  at  a  rectangular  aperture 
y  -  (sin  x/x) 2 


Within  each  rectangle  formed  by  nairs  of 


consecutive  dark  lines  of  the  Fourier  transform  of  the 
aperture  the  intensity  rises  to  a  maximum;  all  their  maxi- 
mums  are,  however,  only  a  small  fraction  of  the  central 
maximum  of  the  Fourier  transform  and  decrease  rapidly  with 
increasing  distance  from  the  center  (see  Fig.  8).  The 
larger  the  rectangular  opening,  the  smaller  the  effective 
size  of  the  diffraction  pattern  of  the  Fourier  transform. 


TABLE  1 

2 

The  First  Five  Maxima  Of  The  Function  y  *  (sin  x/x) 


X 

- 2 - 

y  =  (sin  x/x) 

0 

1 

1.430 

*  4.493 

0. 04718 

2.459 

-  7.725 

0. 01694 

3.470 

=  10.90 

0. 00834 

4.479 

=  14.07 

0. 00503 

Figure  8.  Fraunhofer  diffraction  pattern  of  a  rectangular 
aperture 

b.  Circular  Aperture 

The  Fraunhofer  diffraction  pattern  for  a  cir¬ 
cular  aperture  will  be  investigated.  Using  polar  coordinates 
let  (  p,  6)  be  the  polar  coordinates  of  a  particular  point  in 
the  aperture: 

pcot9<m§,  prinO  *»  r»; 

(36) 


and  let  (  o),\p  )  be  the  coordinates  of  a  point  P  in  the 
diffraction  pattern,  referred  to  the  geometrical  imaqe  of 
the  source: 


w  co tya;,  v  (in  f  ■»  q. 


(37) 
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2  2 

From  the  definition  of  d  and  q,  it  follows  that  w  «  p  +  q 
is  the  sine  of  the  angle  which  the  direction  (o,q)  makes 
with  the  central  direction  p  =  q  *  0.  The  diffraction  in¬ 
tegral  of  equation  (32)  now  becomes,  if  a  is  the  radius  of 
the  circular  aperture, 

U(P)  =  cj^'*-** *»<•-''> pdpdd.  ( 38 ; 

Now  using  the  integral  representation  of  the  Bessel 
function  J  (z) : 


«■  *  e*“  da 


-/.(*>• 


(39) 


Equation  (38)  therefore  reduces  to 

U{P)  -  2ircj*Jo(kp">)pdP- 


(40) 


Also,  the  recurrence  relationship  is 

5  <*"*‘-W*)>  -  **" '«(*)• 


(41) 


giving  for  n  ■  0,  on  integration 


iw*  -  xj^x)- 

'o 


(42) 
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From  equations  (40)  and  (42)  it  follows  that 


(43) 


where  D  *  ir 


Hence  the  intensity  is  given  by 


I(P) 


(44) 


where 

In  =  (PD*  =  ED/ A*. 

Equation  (44)  is  the  formula  derived  by  Airy  in  somewhat 
different  form.  /~13 J 

The  intensity  distribution  in  the  neighborhood 

of  the>  geometrical  image  is  characterized  by  the  function 
2 

Y  38  (2J^(x)/x)  shown  in  Fig.  9a.  This  function  y  has  its 

principle  maximum  y  =  1  at  x  *  0,  and  with  increasing  x,  it 

oscillates  with  gradually  diminishing  amplitude  in  a  similar 

2 

way  to  the  function  (sin  x/x)  which  was  discussed  previously. 
The  intensity  is  zero  (minimum)  for  values  of  x  given  by 
Jx(x)  -  0. 

The  minima  are  no  longer  equidistant  (see  Table  2) . 
The  positions  of  the  secondary  maxima  are  given  by  the  values 
of  x  that  satisfy  the  equation. 
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(a)  The  function 
2J1  (x)  /x)  2 


(b)  Airy  disk 


Figure  9. 


Fraunhofer  diffraction  pattern  of  a  circular 
aperture.  From  Fowles  /  22  7 


TABLE  2 

The  First  Few  Maxima  And  Minima  Of  The  Function 
y  »  ( 2J1  (x ) /x) 2 


x  (2J^  (x)/x)2 


0 

1 

Ma: 

1.220 

m 

3.833 

0 

Mil 

1.635 

m 

5  .136 

0 .0175 

Ma; 

2.233 

m 

7  .016 

0 

Mil 

2  .679 

m 

8  .417 

0 .0042 

Ma: 

3.238 

m 

10  .174 

0 

Mil 

3.699 

m 

11.620 

0 .0016 

Ma: 

or  using  the  equation  [ar"JB(*)i  =  —  (45) 

dx  |_  i 

by  the  roots  of  equations  (x)  *  0  with  increasing 
the  separation  between  two  successive  minima  or  maxima 
approach  the  value  of  it  as  in  the  case  of  the  rectangular 
aperture. 

The  results  show  that  the  pattern  consists  of 
a  bright  disk,  centered  on  the  geometrical  image  at 
p  =*  q  =*  0  of  the  source  and  surrounded  by  concentric 
bright  and  dark  rings  (see  Pigs.  9b  and  10) . 


Figure  10.  Fraunhofer  diffraction  oattern  of  a  circular 
aperture  (the  Airy  pattern)  6  mm  in  diameter, 
magnification  50X,  mercury  yellow  light  X  * 
5790A. 
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The  intensity  of  the  bright  rings  decreases 
rapidly  with  increasing  radius;  only  the  first  or  second 


rings  are  bright  enough  to  be  visible  to  the  naked  eye. 

The  bright  central  area  is  known  as  the  Airy  disk  (see  Fig. 
10) .  The  Airy  disk  extends  to  the  first  dark  ring  whose 
size  is  given  by  the  first  zero  of  the  Bessel  function, 
namely  x  =  3.833.  The  angular  radius  of  the  first  dark 
ring  is  thus  given  by 

sine  *  a  9  (46 ) 

which  is  valid  for  small  angles,  0  .  Here  a  is  the  radius 
of  the  aperture. 

The  image  of  a  distant  point  source  formed  at  the 
focal  plane  of  a  camera  line  is  actually  a  Fraunhofer  dif¬ 
fraction  pattern  for  which  the  aperture  of  the  lens  is  the 
opening.  Thus  the  image  of  a  composite  source  is  a  super¬ 
position  of  many  Airy  disxs.  The  resolution  of  detail  in  the 
image,  therefore,  depends  on  the  size  of  the  individual  Airy 
disks.  If  D  is  the  diameter  of  the  lens,  then  the  angular 
radius  of  the  Airy  disk  is  approximately  1.221/d.  This  is 
also  the  approximate  minimum  angular  separation  between  two 
equal  point  sources  such  that  they  are  just  barely  able  to 
be  resolved,  because  at  this  angular  separation,  the  central 
maximum  of  the  image  of  one  source  falls  on  the  first  minimum 
of  the  other  (see  Fig.  11) .  This  condition  for  optical  reso¬ 
lution  is  the  Rayleigh  criterion.  The  Rayleigh  criterion 
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does  not  always  give  an  absolute  minimum  separation  for 
resolution,  but  it  is  almost  always  used  because  of  its 
simplicity.  /T3 ,  17,  22,  527 

Saddle  point 


Figure  11.  The  Rayleigh  criterion. 

D.  FOURIER  TRANSFORMS 
1.  Convolutions 

Most  considerations  of  diffraction  involve  the  use 
of  the  Fourier  transform  in  one  form  or  another.  One  of  the 
most  important  properties  of  the  Fourier  transform  involves 
the  concept  of  a  convolution,  or  convolution  integral; 
and  for  both  this  and  the  Fourier  transform,  it  is  convenient 
to  use  delta- functions. 

Consider  a  special  delta  function,  the  Dirac  delta 
function.  A  Dirac  delta  function  at  x  ■  a  is  defined  by 

(0  for  x*a 
"  foi  **« 

and  (47) 

as 

J  &(x-a)dx  •  1  , 

—  m 

The  delta  function  at  x  ■  0,  5  (x) ,  can  be  con¬ 

sidered  as  the  limit  of  a  set  of  real  continuous  functions. 
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such  as  Gaussians: 


swMr.[^;ex|,(-vi]-  ms) 

As  a  tends  to  infinity,  the  Gaussian  function  has 
a  maximum  value  tending  to  infinity  and  a  half-width  (1/a) 
tending  to  zero,  but  the  integral  over  the  function  is  always 
unity.  The  delta  function  can  be  used  as  a  convenient  no¬ 
tation  for  any  function  of  integral  unity  in  the  form  of  a 
sharp  peak  having  a  width  so  small  it  is  not  experimentally 
significant. 

The  definition  of  a  delta  function  is 


5(x)*  J  exp{2ffurj'}<l>' ,  (49) 

—  m 

*  In  one  dimension,  the  convolution  integral  (or  con¬ 
volutions  product)  of  two  functions  f (x)  and  g(x)  is  defined 
as 

•• 

C(x)*f(x)>f(x)s  f  f(X)g(x-X)dX .  (50) 

By  a  simple  change  of  variable. 


/(*)•*<*)-  /  g(XMx-X)dX*g(x)>f(x). 


(51) 


The  identity  operation  is  the  convolution  with  the 


Dirac  delta  function: 


f{x)  *  5(jc)  =/(.t)  , 

/(x)  *  $(x  -  a)  =  /(x  -  a)  . 


(51a) 


For  a  general  intensity  distribution,  I (x) ,  each 
sharp  spectral  line  or  each  part  of  a  broader  spectral  line 
will  be  "spread  out"  by  the  "spread-function"  g(x),  so  that 
the  recorded  intensity,  I0BS  be  ^ess  sharply  peaked  or 

less  well  resolved  than  the  original  spectrum  (see  Fig.  12). 


I(x)jr<x-X) 


C(x)  =  I(x)*g(x> 


Figure  12.  A  convolution  operation.  An  intensity  function 
I(x)  representing  a  spectral  line  is  multiplied 
by  the  transmission  function  g(x-X)  of  a  slit, 
centered  at  X*x;  the  product  of  these  two  functions 
is  integrated  to  give  the  measured  intensity  I0BS (x) . 

Likewise,  the  blurring  of  an  image  due  to  the  im¬ 
perfections  of  a  camera  lens  may  be  described  in  terms  of 
convolution  of  the  ideally  perfect  image  intensity  with 
some  two  dimensional  function  g(x,y).  For  a  point  source 
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of  light,  the  ideal  image  would  be  a  two  dimensional 
function  6(x,y).  The  spreading  gives 


A>bi  *t(x,y)  •  5(*.j0 s  g(x<y)  ■ 


(52) 


For  a  general  object,  consisting  of  a  large  number 
of  independently  emitting  point  sources  and  having  an  ideal 
image,  I0(x,y), 

Iqbs  =  9(x»y)  *  Io(x'y)  (53) 


Thus  each  point  of  the  original  intensity  distri¬ 
bution  is  spread  into  a  disk  of  intensity  of  light  and  the 
overlapping  of  these  disks  gives  a  blurring  and  loss  of 
resolution  of  the  image. 

The  Fourier  transform  of  a  one-dimensional  function 
f (x)  is  defined  as 


7  [fix)]  s  F(u)  *  f  Ax)  exp  (2rriux}dx . 


(54) 


The  inverse  transform. 


is  defined  so  that 


/(x)  -  7 _1  [?{/(*»]  a  J  F{u)exp{-2mux}4u.  (55) 
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Here  the  convention  of  including  2ir  in  the  exponent 
is  followed.  This  convention  is  commonly  used  in  considera¬ 


tions  of  diffraction  in  order  to  avoid  the  necessity  of 
adding  a  constant  multiplier  2tt  in  equations  (54)  or  (55J./X77 
2.  Diffraction  Analysis 

The  amplitude  of  scattering  from  an  object  in  the 
Fraunhofer  diffraction  approximation,  whether  derived  from 
the  Kirchhoff  formulation  or  from  scattering  theory,  is  de¬ 
scribed  by  a  Fourier  transform  integral.  Thus  the  diffraction 
amplitude  may  be  described  in  terms  of  a  distribution  in 
Fourier  transform  space,  often  referred  to  as  "reciprocal 
space . " 

'f  The  individual  properties  of  Fourier  transforms  will 
not  be  expounded  on  here,  as  they  are  well  documented  in 
many  references  readily  available.  The  remainder  of  this 
section  will  demonstrate  the  application  of  Fourier  trans¬ 
forms  and  their  relationship  to  diffraction. 

The  amplitude  distribution  of  a  very  small  source  or 
the  transmission  through  a  vary  small  aperture  (or  slit)  in 
one  dimension  may  be  described  as  6  (x)  or  by  <5(x-a)  when  it 
is  not  at  the  origin.  The  Fourier  transform  used  to  derive 
the  Fraunhofer  diffraction  pattern  is 


t  ?$(*)»  I, 

76(x-a)  *  exp{2niua}  . 


(56) 


To  illustrate  the  above,  the  following  integral  is  written: 
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(57) 


J  6(x -a)  exp{2jriux}dx  . 

—  a* 

The  integrand  is  zero  except  for  x  *  a.  Hence  the  integral 
may  be  written 


exp{2ri ua}  f  S(x-a)dx  <*  exp(2triua}  . 

The  amplitude  of  the  diffraction  pattern  will  be 

proportional  to  F  (u)  -  F  5 (x-a)  where  u  =  1/A  .  The 

2 

intensity  observed  will  be  proportional  to  F(u)  *  1. 

The  Fourier  transform  of  a  plane  wave  with  respect 
to  t  only,  gives 


7,  exp{2jri(>V-*/*i)}  a«(‘,+>'i)  exp{-2iri x/At>  ,  (59) 

i.e.,  the  delta  function  in  frequency  and  a  plane  wave  in 
real  space.  Fourier  transform  with  respect  to  x  gives 


9X  exp{27ri(i>j  f-x/Aj)}  ■  exp{2nv,  r} 5(«-  1/Xj) ,  ( 6  ®  ^ 

i.e.,  a  sinuisoidal  variation  in  time  and  a  delta  function 
in  reciprocal  space.  Transforming  with  respect  to  both 
variables  gives 
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9x  t  exp  {2rti(vl  t  -  x /A  i ) }  =  S(»'+«'i ) •  5(u  - J  /Xt )  , 


(61) 


which  is  a  delta  function  in  both  frequency  and  reciprocal 
space . 

The  Fourier  transform  of  all  possible  waves  which 
can  exist  in  a  given  medium  with  respect  to  both  x  and  t 
is  a  set  of  points  defining  the  relationship  between  the 
frequency  v  and  (or  between  the  angular  frequency 

u  and  K)  which  is  known  as  the  dispersion  relation  for  that 
type  of  wave  and  the  particular  medium. 

For  the  translation  of  an  object,  the  Fourier  trans¬ 
form  is  written, 

<7 fix -a)  =■  71 /(*)  *  6(x-tf)l  -  exp  {2iriua}  ( 6  2 ) 

Thus  the  translation  of  the  object  in  real  space  has  the 

effect  of  multiolying  the  amplitude  in  reciprocal  soace  by 

a  complex  exponential.  The  intensity  distribution  of  the 

2 

Fraunhofer  diffraction  pattern  is  given  by  F(a)  ,  which 
is  independent  of  translation. 

a.  Rectanqular  And  Circular  Apertures 

The  intensity  distribution  for  both  the  rec¬ 
tangular  and  circular  apertures  is  important  when  studying 
Fraunhofer  diffraction.  The  transmission  function  of  a 
rectangular  aperture  in  two-dimensional  form  may  be  written 
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f(x.y) 


i 

o 


if  |x|<a/2  and  \y\<bil. 


(63) 


Then 


j/2  i/2 

F{u,v)=  J  exp{2ffiux}dx  J  exp{2ffioy}dj?'  , sin(irgu)  sin (nbv) 
-i/2  -i/2  nau  nbv  ' 


(64) 


so  that/  for  diffraction  from  a  rectangular  aperture,  the 
intensity  distribution  is 


l(u,v)-a2b2 


sin2  (tom)  sin2(ffit>) 
(nau)2  (nbv)1 


(65) 


,  2  2 
The  maximum  intensity  atus'v^Oisab. 

Zeros  occur  at  intervals  of  a”^  along  the  u  direction, 

which  is  parallel  to  the  x  axis  and  at  b*1  intervals  along 

the  v direction  (see  Fig.  13). 


b  .. 

a 

q(x,y) 


■e* 


0  llll  0 

F(u,v) 


Figure  13.  Diagram  suggesting  the  form  of  the  Fourier 
transform  of  a  rectangular  aperture.  From 
Cowley  /""17  7 
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Thus  the  intensity  function  falls  off  through  decreasing 
oscillations  in  each  direction.  The  dimensions  of  the  dis¬ 
tribution  in  reciprocal  space  are  inversely  proportional 
to  the  dimensions  in  real  space. 

The  transmission  function  for  a  circular 
aperture  is  written 

1  if  (j^+y2)1'1  <a/2  , 


f(x,y) 3 1 


(66) 


The  Fourier  transform  is  best  carried  out  by 
writing  the  Fourier  integral  in  plane  polar  coordinates. 

The  result  is: 

w  ,  AO"")  (67) 

Ftu)  3  a  - -  '  ' 

nau 

where  u  is  the  radial  coordinate  and  J^(x)  is  the  first  order 
Bessel  function.  The  function  J^(x)/x  is  similar  to  the  form 
(sin  x)/x  but  has  a  somewhat  broader  central  maximum  with 
the  first  zero  at  1.22a"1  instead  of  a”1. 

b.  Single  And  Double  Exposure  Specklegrams 

Before  calculating  the  diffraction  pattern  from 
a  double -exposure  specklegram,  first  consider  a  single-ex¬ 
posure  photographic  transparency  of  an  object  illuminated 
by  laser  light.  The  recorded  image  is  speckled.  If  a  small 
diameter  laser  beam  is  passed  through  a  portion  of  the 
transparency,  the  aperture  w  will  be  the  amplitude  trans¬ 
mission  function  F.(x,y)  inside  the  illuminated  area.  Thus 


cF(R(x,y) 


(68) 


u 


and 


io«  p[R(x,y)]  *  F  ^R(x,y)]  (69) 

Now  for  a  double  exposure  speckle  photographic 
transparency  where  the  object  has  suffered  an  in-plane 
displacement  ax  in  the  x  direction  between  exposures,  the 
amplitude  transmission  function  is 


R2(x,y)  *  R (x  -  Ax/2/y)  +  p(*  +  Ax/2, y)  (70) 

A  series  expansion  of  equation  (70)  yields 

R,(x,y)  -  2  feu.y)  *  (il)2  +  (0*) 

1  L  21  '  2'  41  V  2J 


+  .  .  . 


(71) 


where  the  prime  indicates  a  differential  with  respect  to 
the  variable  x.  The  linear  prooerties1  of  a  Fourier  transform 
provide  that 


F 

[gn  (x)]  - 

(iw  )n  F  [g(x)J 

(72) 

where 

gn(x)  is 

the  nfc^  derivative 

of  the 

function  g(x). 

Thus 

U 

”  c  f[r2 

(x,y)j 

U 

»  c  fJr(x,v)J2 

i  1  Ax 

1  "  jr  w,-^. 

L 

+ 

i 

CM 

Ax  4  1 

w,i=2  +... 

U 

*  c  F  [R(x,y)J 

2  cos  (w, 

) 

(73) 

(73) 


Using  equations  (72) and  (73)  and  w  *  K  sine 
U  =  UQ  2  cos(  ^  K  sin6  )  (74) 

I  *  IQ  +  cos  (KA  x  sine  )j  (75) 

Thus  the  intensity  distribution  for  a  specklegram 
is  the  diffraction  halo  from  a  single -exposure  transparency 
modulated  by  Young's  fringes.  /~3 , A7 

3.  Two  Dimensional  Fourier  Transforms 

Two-dimensional  Fourier  transforms  are  fundamental 
to  all  image  formation.  The  Fraunhofer  diffraction  pattern 
of  an  object  is  the  Fourier  transform  of  that  object, 
object  of  low  spatial  frequency,  i.e.,  an  object  that  has 
large  details,  and  wide  separation  between  details,  will  give 
a  Fourier  transform  of  small  elements  that  are  close  together. 
An  object  of  high  spatial  frequency  will  give  maxima  that  are 
wide  apart. 

a.  Thickness  Function 

To  find  the  thickness  A(x,y)  of  a  lens  at  any 
point  (x,y)  we  split  the  lens  into  two  parts,  as  shown  in 
Fig.  14,  and  write  the  total  thickness  function  as  the  sum 
of  the  two  individual  thickness  functions, 

A  (x,y)  *  Aj(x,y)  +  A2  (x,y)  (76) 

Referring  to  the  geometry  shown  in  Figure  (14),  the  thick¬ 
ness  A^(x,y)  is  given  by 
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&\(x,y) 


-  («,  -  VRx i  -  *»  -  y'-) 


(77) 


Similarly, 


A»(x,#)  *  A**  —  (—Rt  —  ^ Rt*  ~  x*  ~  y*) 

-  Ao*  +  *«(i  -  V1  £r) 


(78) 


where  we  have  factored  the  positive  number  out  of  the 
square  root.  Combining  equations  (76),  (77),  and  (78),  the 
total  thickness  function  is  seen  to  be 


(79) 


Figure  14.  Calculation  of  the  thickness  function.  From 
Goodman  /"* 2  5_7 
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b.  Paraxial  Aporoximation 

The  expression  for  the  thickness  function  can 
be  substantially  simolified  if  attention  is  restricted  to 
portions  of  the  wave  front  that  lie  near  the  lens  axis,  or 
equivalently,  if  only  paraxial  rays  are  considered.  Thus  we 
consider  only  values  of  x  and  y  sufficiently  small  to  allow 
the  following  approximations  to  be  accurate: 


l  _ *L+y? a I  _ x'  +  yt 


R  i* 


X1  +  y* 

1  -  -  =.  1  - 
Rtl 


2Ri* 
x'  +  y* 

2Ri* 


(80) 


The  resulting  phase  transformation  will,  of  course,  represent 
the  lens  accurately  over  only  a  limited  area,  but  this 
limitation  is  no  more  restrictive  than  the  usual  paraxial 
approximation  of  geometric  optics.  Note  that  equation  (80) 
amounts  to  approximation  of  the  spherical  surfaces  of  the 
lens  by  parabolic  surfaces.  With  the  help  of  these  approxi¬ 
mations,  the  thickness  function  becomes 

*«,-*-£TJi(s;-s).  <81) 


c.  Phase  transformation 

The  paraxial  approximation  to  the  lens  trans¬ 
formation  is: 

*»(*.?)  -  exp  (JknAo]  exp  £  -jk(n  -  1)  -  t  -L* 
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where  t^(x,y)  is  the  lens  transformation  function.  The 
physical  properties  of  the  lens  (that  is  n,  ,  and  R£)  can 
be  combined  in  a  single  number  f  called  the  focal  length; 
which  is  defined  by 

-f&  (n  -  1)  (jf,  "  r-)  (82) 

Thus  the  phase  transformation  may  be  rewritten 

,  $x,y)  -  exp  L?fc«A.]  exp  £ -j  ~  (**  +  Vs)]  (83) 

This  equation  will  serve  as  our  basic  representation  of  the 
effects  of  a  lens  on  an  incident  disturbance. 

One  of  the  remarkable  and  useful  properties  of 
a  converging  lens  is  its  inherent  ability  to  perform  two 
dimensional  Fourier  transformations. 

d.  Fourier  Transform  Configurations 

Three  separate  configurations  for  performing 
the  Fourier  transform  operation  are  considered  (see  Fig.  15). 


Object 


(a)  Object  placed  against  the  lens. 

Fioure  15.  Fourier  Transforming  Configurations. 


(b)  Object  placed  in  front  of  lens. 


(c)  Object  placed  behind  lens. 

Figure  15.  (Continued). 

In  the  material  that  follows,  the  illumination 
is  assumed  monochromatic  and  the  distribution  of  light 
amplitude  across  the  back  focal  plane  of  the  lens  is  of 
concern.  The  front  and  back  focal  planes  of  a  lens  are 
defined  as  the  planes  normal  to  the  lens  axis  situated  at 
a  distance  f  in  front  of  and  in  back  of  the  lens  re¬ 
spectively.  The  front  of  the  lens  is  the  side  of  the  lens 
the  incoming  beam  strikes  first. 

(1)  Object  Placed  Against  The  Lens.  For  the  first 
case  where  the  object  is  placed  against  the  lens,  let  a  plane 
object  with  amplitude  transmittance  tQ(x,y)  be  placed  imme¬ 
diately  in  front  of  a  converging  lens  of  focal  length  f,  as 
shown  in  Fig.  15a.  The  object  is  assumed  to  be  uniformly  il¬ 


luminated  by  a  normally  incident,  monochromatic  plane  wave  of 


amplitude  A,  in  which  case  the  disturbance  incident  on  the 
lens  is 


Urfoy)  -  AK(x,y) 


The  finite  extent  of  the  lens  aperture  can  be  accounted  for 
by  associating  with  the  lens  a  pupil  function  P(x,y)  defined 


P(x,v)  =  I  1  inside  the  lens  aperture 

1 0  otherwise 


Thus  the  amplitude  distribution  behind  the  lens  becomes, 
using  equation  (83), 


U'»(*,V)  m  Ui (x,y)P(x,y)  exp  ^  (**  +  y*)j 


The  constant  phase  delay  associated  with  the  lens  transfor¬ 
mation  has  been  omitted  since  it  does  not  affect  the  result 
in  any  significant  way. 

To  find  the  distribution  Uf(Xfy^)  of  field 
amplitude  across  the  back  focal  plane  of  the  lens,  the  Fresnel 
diffraction  formula  is  applied. 

Thus,  putting  z  «  f. 


U/(x/,y,) 


“p  [;  (Z/’  +  y/)]  ft  , 

exp  (x*  +  y!)j  exp  ^  (xx,  +  yy/)]  dx  dy 


where  a  constant  phase  factor  has  been  dropped.  Substituting 
equation  (86)  in  equation  (87),  the  quadratic  phase  factors 
within  the  integrand  are  seen  to  cancel,  leaving 


U/(*/.y/) 


/ f  Ui(z,y)P(x,y)  exp  J^-y  ^  {xx,  +  yyf)  j  dx  dy 


(88) 


Thus  the  field  distribution  is  proportional  to  the  two- 
dimensional  Fourier  transform  of  that  portion  of  the  incident 
field  subtended  by  the  lens  aperture.  VThen  the  physical 
extent  of  the  object  is  smaller  than  the  lens  aperture,  the 
factor  P(x,y)  may  be  neglected,  yielding 

U t(x„yt)  -  il— I V'  ^  //  6Xp  [~J1j  (zx'  +  W)]<£*dy  <"> 

Evidently  the  amplitude  and  phase  of  the  light  at  coordinates 
(Xf,yf)  are  influenced  by  the  amplitude  and  phase  of  the 
object  Fourier  component  at  spatial  frequencies 
(fx  ■  xf/*f  fy  =*  y^/Af ) .  Measurement  of  the  intensity 
distribution  yields  knowledge  of  the  power  spectrum  of  the 
object;  the  phase  distribution  is  of  no  consequence  in  such 
a  measurement .  Thus 


(90) 


(2)  Object  In  Front  Of  Lens.  Consider  next 
the  more  general  geometry  of  Fig.  15b.  The  object,  located 
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a  distance  dQ  in  front  of  the  lens,  is  illuminated  by  a 

normally  incident  plane  wave  of  amplitude  A.  The  amplitude 

transmittance  of  the  object  is  again  represented  by  tQ.  In 

addition,  let  FQ(£^,f  )  represent  the  Fourier  spectrum  of 

the  light  transmitted  by  the  object,  and  F  (f  ,f  )  the  Fourier 

A  X  Y 

spectrum  of  the  light  incident  on  the  lens;  that  is. 


Ufxjr)  -  ff{.4t.landF1(/r,/r)  -  5IU‘>  {9] 

fiifx.fr)  -  f. ifx.fr)  exp  f  —jr\d.(fx*  +/r*)  1 

For  the  moment,  the  finite  extent  of  the 
lens  aperture  will  be  neglected.  Thus,  letting  P  *  1, 
equation  (88)  can  be  rewritten 


U/(*/,y/) 


[;  |  (V +  »,•)]  .  . 

jV  '  W  W 


Substitution  of  eouation  (91)  in  ecruation  (92),  gives 

... ... 


U/(x/,y/) 


-  (**  yj\ 

f-  W  x7j 


U/(x/,y/) 


A  cxp  [;  2f  Q  +  v,t^\  //  *•(*«»•)  oxp  [  -i  ^  (*.*/  +  j t>y/)  ]  (93) 


Thus  the  amplitude  and  phase  of  the  light  at  coordinates 
(Xj,Yj)  are  again  related  to  the  amplitude  and  phase  of  the 
object  spectrum  at  frequencies  (x^  /Af,  y^/Af  ).  Note  that 
a  phase  factor  again  precedes  the  transform  integral,  but 
that  it  vanishes  for  the  very  special  case  d  ■  f. 

(3)  Object  Placed  Behind  The  Lens.  Finally,  we 
consider  the  case  of  an  object  that  is  placed  behind  the 
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lens,  as  shown  in  Fig.  15c.  The  object  again  has  trans¬ 
mittance  tQ,  but  is  now  located  a  distance  d  from  the  focal 
plane  of  the  lens.  Let  the  lens  be  illuminated  by  a  nor¬ 
mally  incident  plane  wave  of  amplitude  A.  Then  incident  on 
the  object  is  a  spherical  wave  converging  toward  the  back 
focal  point  of  the  lens. 

In  the  geometrical-optics  approximation,  the 
amplitude  of  the  spherical  wave  impinging  on  the  object  is 
Af/d;  the  particular  region  of  the  object  that  is  illuminated 
is  determined  by  the  intersection  of  the  converging  cone  of 
rays  with  the  object  plane.  If  the  lens  is  circular  of 
diameter  l,  then  a  circular  region  of  diameter  £  d/f  is  illumi¬ 
nated  in  the  object  space.  The  finite  extent  of  this 
illuminating  spot  can  be  represented  mathematically  by  pro¬ 
jecting  the  pupil  function  of  the  lens  along  the  cone  of 
rays  onto  the  object  yielding  an  effective  oupil  function 
P  |xQ (f/d)  ,yQ (f/djj  in  the  object  plane.  Note  that  the  object 
transmittance  tQ  will  have  a  finite  aperture;  the  effective 
aperture  in  the  object  space  is  thus  determined  by  the  inter¬ 
section  of  the  true  object  aperture  with  the  projected  pupil 
function.  If  the  object  aperture  is  fully  illuminated,  then 
the  projected  pupil  function  may  be  ignored. 

Using  a  paraxial  approximation  to  the  spheri¬ 
cal  wave  that  illuminates  the  object,  the  field  amplitude 
transmitted  by  the  object  may  be  written 
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1 


U.(x.,y.) 


(94) 


Assuming  Fresnel  diffraction  from  the  object  plane  to  the 
focal  olane 


(95) 


Thus,  up  to  a  quadratic  phase  factor,  the  focal-plane 
amplitude  distribution  is  the  Fourier  transform  of  that 
portion  of  the  object  subtended  by  the  projected  lens  aper¬ 
ture. 

The  result  represented  by  equation  (95)  is 
essentially  the  same  result  obtained  previously  when  the 
object  was  placed  directly  against  the  lens  itself.  How¬ 
ever,  an  extra  flexibility  has  been  obtained  in  the  present 
configuration;  namely,  the  scale  of  the  Fourier  transform 
is  under  the  control  of  the  experimenter.  By  increasing  d, 
the  spatial  size  of  the  transform  is  increased,  at  least  un¬ 
til  the  transparency  is  directly  against  the  lens  (i.e.,  d  *  f). 
By  decreasing  d,  the  size  of  the  transform  is  made  smaller.  /75_7 

E.  INTERFERENCE 

The  phenomenon  of  interference  is  central  to  the  subject 
matter  presented  within.  To  initiate  consideration  of  inter¬ 
ference,  two  different  light  waves,  and  E2,  of  the  same 
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frequency  are  superimposed.  Since  E  *  +  E2,  the  irradiance 

will  be: 


I  -  <E2>  -  <E2>  +  <Et>  +  2<E,.Ei>.  ( 96 ) 

where  <>  symbolizes  the  time  average  or  mean.  The  pro¬ 
portionality  constant  ev  of  equation  (5)  will  be  neglected 
for  the  remainder  of  the  discussion.  From  the  previously 
stated  assumption,  both  waves  are  linearly  polarized  in  the 
same  direction.  Then 


Ei  —  A  |  cos  (<o/ — k|  -r) 


(97) 


E2  m  A  j  cos  («/  —  k2*r + <£), 


(98) 


where  0  is  a  constant  relative  phase  between  the  two  waves, 
and  and  K2  are  the  wave  propagation  vectors  for  and 
E2  respectively.  Combining  equations  (96),  (97),  and  (98) 
and  carrying  out  the  averaging,  then 


/"/,  +  /,  +  2  V/ITj  cos  8 


(99) 


where 


where  /,"‘Aft/2mA J, 


and 


6  -  kj»r  -  k,f  -  ^ 


(100) 


(101) 
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where  a  is  the  phase  difference  between  the  two  waves  at  any 
location.  The  irradiance  varies  from  a  minimum  value 
IMin  "  +  *2  “  2(iiI2^'S  at  P°ints  where  <5  =  (2N  +  1) 

j- 

to  a  maximum  value  IMAX  *  1^  +  I2  +  2^I1I2^  at  P°^nts 

where  <5  *  2Nir  ,  where  N  is  an  integer.  The  irradiance 

pattern  will  be  alternating  light  and  dark  fringes. 

1.  Interference  From  Two  Point  Sources  Emitting  Spherical 
Waves 

A  specific  important  example  of  interference  is  that 
of  two  spherical  waves  emanating  from  two  point  sources  of 
light  (see  Fig.  16).  Assume  that  light  is  radiated  from 
points  b  and  c  in  phase,  that  is  6  3  0.  The  irradiance  at 
any  point  Pn  in  space  is  given  by. equation  (99)  with  the  phase 
difference  being  f. 


Figure  16.  Interference  of  light  emitted  by  two  point 

sources,  and  S2.  Modified  from  Vest  /” 5 3_7 

If  l  is  the  perpendicular  distance  from  the  axis  formed  by 
S^S2  to  the  point  Pn  and  D  ■  d/2  is  the  distance  S^S2,  then 


=*VA  +  (d  -  x)  (1C 

r2  -V^2  ♦  (d  +  x)2  (ic 

The  phase  at  due  to  and  S2  respectively  is 

Kx1  »  2-nx^X  (10 

Kr2  »  2irr2/A  (10 

where  K  =  2ir/A 


The  phase  difference  generated  at  the  point  P  b\ 

n 


and  S2  is  then 


6  -  K(rx  -  r2) 


Substituting  the  values  for  r,  and  r~ 


using  the  power  identity  of 


(1  +  x)n  -  1  +  nx  +  +  ...  (109) 

and  neglecting  higher  order  terms 

6  »  2Kxd/ 1  (110) 

Then  the  locus  of  points  forming  a  surface  of  maximum 

irradiance  is  determined  bv  setting  ■  2wN 

where 

r.  -  r2  *  m  n  X  ;  N  -  0,  1,  2  ...  (Ill) 

This  is  the  equation  of  a  family  of  hyperboloids 
of  revolution  about  the  axis  S^sJ  connecting  the  two  point 
sources. 

Two  light  waves  which  are  capable  of  interfering 
are  said  to  be  coherent.  Because  of  the  coherence  require¬ 
ments,  most  interference  experiments  are  conducted  using 
two  images  of  the  same  physical  source. 

2.  Young's  Fringes 

A  simple  division  of  wave  front  interferometer  is 
used  to  produce  Young's  fringes.  This  is  shown  in  Fig.  17  . 
It  is  simply  an  opaque  screen  in  which  two  holes  separated 
by  a  distance  0  have  been  cut.  This  screen  is  illuminated  by 
a  point  source  located  a  distance  £g  from  the  screen  and  a 
distance  y  above  the  center  line  axis  of  symmetry  of  the 

9 

screen.  The  light  diffracted  by  the  two  holes  forms  an 
interference  pattern  which  can  be  observed  on  a  screen  placed 
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some  distance  £  away.  In  practice,  y  ,  y,  and  D  are  much 

O  5 

smaller  than  ' i  and  £  . 

s  o 


Figure  17.  Interferometer  to  form  Young's  fringes. 

The  irr&diance  of  each  of  the  light  waves  at  y 
therefore  will  be  nearly  equal  1^  =  I2  *  IQ  so  equation  (99) 
becomes 


I  *  21  (1  +  cos  a) 

o 

I  -  4Iq  cos2  (  a/2) 


(112) 


Here  the  phase  difference  is  6  **  KA£,  where  A£ 
is  the  difference  in  the  distance  the  light  travels  from 
the  source  S  to  the  observation  screen  at  y  on  the  two 
paths. 


A£ 


*  <d  -  v 2] 11  * 

[tsJ  ♦  M  *  y,)2j  15  - 


[t02  *  M  -  y)2]  11  - 
[t02  ♦  (d  +  y)  2j  h 


(113) 
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Dividing  by  lg  and  lQ  to  set  the  equation  up  to  utilize 
the  identity  of  equation  (109)  and  disregarding  second 
order  and  higher  terms  gives 

2-i ) 


At 


V 

s 


At 


'r 

1  +  i 

(d-ys)2 

1 

1  «. 

A  +  7 

1 

t  2 
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Li 
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7 

t  2 

o  J 

X 
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hi 


t>y  Dy 

s  o 


(114) 


Knowing  that  6  »  KAt  gives 


6  »  -K 


Dys  Dy 

T —  +  T~ 
s  o 


;] 


(115) 


Where  again  K  *  2*  /X  is  the  wave  number.  The 
irradiance  is  then  equal  to 


4I0  cos' 


0) 


(  116) 


Substituting  for  5  into  equation  ( 116)  gives 


(117) 


I  =  41  cos' 
o 


jl°  [JL  +  IS. 

\  U.  i. 


which  represents  parallel  fringe  spacing  of  XI ^D. 

3.  Scalar  Wave  Theory 

Further  consideration  of  the  subject  of  diffraction, 
shows  the  representation  of  light  waves  may  be  modified. 

Seeing  that  the  assumption  of  unidirectional  polarized  light 
has  been  made,  only  one  vector  component  of  E  need  be  studied. 
In  other  words,  E  can  be  treated  as  a  scalar  instead  of  a 
vector.  E  is  then  replaced  by  the  scalar  quantity  U(x,y,z,t), 
which  is  referred  to  as  the  optical  disturbance.  For  example, 
a  monochromatic  plane  wave  is  denoted  by 


U (x,y,z,t)  »  a(x,y,z)  cos(wt  -  K*r) 


(118) 


Similarly,  the  expression  for  any  momochromatic  wave  can  be 


written  in  the  form 


U(x,y,z,t)  -  a(x,y,z)  cos  |wt  -  0  (x,y,z)j 


(119) 


where  a(x,y,z)  is  the  real  amplitude  of  the  light  wave,  and 
9  (x,y,z)  is  its  phase.  Of  course,  the  plane  wave  is  a 
special  case  for  which  0  (x,y,z)  ■  K*r  is  constant.  Equation 
(119)  also  can  be  written  as 


where 


U(x,y,z,t)  ■  Re  iu(x,y,z)  exp  (iwt)l 


lf(x,y ,z)  ■  a (x,y ,z)  exp 


j*-i  9(x,y,z) 


(120  ) 


h 


1 


is  called  the  complex  amplitude  of  the  light.  Because  the 
temporal  frequency  of  light  is  so  high,  there  is  no  need  to 
consider  the  time  dependence  ext(-iwt)  explicitly.  The  com¬ 
plex  amplitude  U(x,y,z)  contains  all  the  information  about 
the  spatial  structure  of  liqht  waves  essential  to  speckle 
photography . 

The  complex  amplitude  of  light  emitted  from  a  rough 
surface  received  at  a  detector  perpendicular  to  that  surface 
is  given  by  equation  (120). 

The  irradiance  of  the  light  at  the  detector  surface 
is  a  more  important  quantity  because  that  is  what  is  actually 
detected.  The  probability  density  function  for  the  irradi¬ 
ance,  Pj  (I),  is 

Pt  (I)  -  ^  exp  (121) 

The  irradiance  of  the  speckle  pattern  there  obeys 
negative  exponential  statistics.  A  measure  of  the  contrast 
of  speckle  is  the  ratio  C  ■  o^/  <I>  ,  where  I  is  the 
standard  deviation  of  the  irradiance,  and  I  is  its  mean 
value.  For  the  distribution  given  by  equation  (121),  the 
speckle  contract  is  unity.  /767 

The  most  probable  irradiance  in  a  speckle  pattern 
is  seen  to  be  zero,  that  is  black  (see  Fig.  18). 

F.  SPECKLE  THEORY 

The  speckle  phenomenon  has  long  been  familiar,  but  only 
the  introduction  of  lasers  has  brought  a  deeper  understanding 
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and  many  new  applications  to  this  field.  /” 2 1_7  Laser 
speckles  appear  whenever  an  optically  rough  surface  is  il¬ 
luminated  with  highly  coherent  light.  The  microscopic, 
random  height  variation  of  roughness  of  the  surface  needs  to 
be  of  the  order  of  the  wave  length  of  light  used  to  illumi¬ 
nate  the  object  surface  or  coarser.  This  encompasses  most 
surfaces  except  those  highly  polished.  The  optical  wave 
originating  from  any  moderately  distant  point  consists  of  a 
series  of  coherent  waves,  each  arising  from  a  different  rough 
element  of  the  surface.  The  optical  path  difference  of  these 
various  waves  may  be  several  wave  lengths.  These  waves 
scattered  from  the  optically  rough  surface  have  not  only  dif¬ 
fering  phases  but  also  differing  wave  amplitudes.  The  inter¬ 
ference  of  the  dephased,  but  coherent,  secondary  waves  results 
in  a  granular  pattern  of  intensity  called  speckle. 

1.  Speckle  Characteristics 

Consider  the  factors  which  quantitatively  determine 
the  characteristics,  size,  shape,  and  intensity  form  of  the 
speckle. 

Figures  (18  a  -f  )  are  photographs  of  a  tyoical 
speckle  pattern  produced  by  scattering  light  from  a  rough 
surface. 

The  physical  explanation  of  speckle  is  that  each  point 
on  the  surface  scatters  the  incoming  laser  light  in  such  a 
manner  that  the  scattered  light  from  one  point  interferes  with 
that  of  another  point.  The  randomness  of  the  speckle  pattern 
itself  is  generated  by  the  surface  roughness  which  generates 
phase  differences  in  the  reflected  waves  upon  scattering. 


Micrographs  of  speckle  for  various  lens  aperture 
settings,  showing  the  speckle  size  dependence  on 
imaging  aperture.  Image  magnification  equaled 
.137.  Microscopic  magnification  300X. 


(c)  Speckle  size  imaged  by  lens  of  f/4 


(d)  Speckle  size  imaged  by  lens  of  f/5,6 
Figure  18.  (Continued) 


79 


(e)  Speckle  size  imaged  by  lens  of  f/8 


(f)  Speckle  size  imaged  by  lens  of  f/16. 


Figure  18.  (Continued) 
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To  determine  the  statistical  characteristics  of  the 


irradiance  in  the  observation  plane  of  the  speckle,  it  is 
convenient  to  consider  an  ensemble  average  at  a  single  detec¬ 
tion  point  rather  than  the  whole  field  of  view.  It  is  assumed 
that  the  incident  light  is  monochromatic  and  linearly  polar¬ 
ized,  and  that  its  polarization  is  not  changed  by  the  scatter¬ 
ing. 

The  most  important  statistical  characteristic  of  laser 
speckle  is  its  size.  From  equation  (117),  it  is  seen  that  any 
two  points  on  the  object  surface  separated  by  a  distance  £ 
give  rise  to  fringes  of  frequency,  f  =  l/X z.  The  average 
fringe  frequency  will  be 


1  7  £max  *  T  <L^> 


(122) 


where  L  is  the  length  of  the  illuminated  area  in  the  y 
direction  and  z  is  the  distance  from  the  illuminated  object 
to  the  image  plane.  Hence  the  irradiance  distribution  across 
a  "typical  speckle"  will  be 

(123) 


I(y) 


1  +  cos  2  tt (Ly/3  Xz) 


where  y  is  the  vertical  axis  value. 

The  width  of  the  speckle  can  be  described  as  the  dis¬ 
tance  between  points  where  I  drops  off  to  one-half  its  maxi¬ 
mum  value.  This  width  is  1.5  (X  z A)?  then  a  typical  speckle 
width  b_  would  be 

3 


b 


s 


-  1.5 


(124  ) 
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To  estimate  the  speckle  size  in  this  case,  treat  the  disk 
enclosed  by  the  pupil  of  the  lens  as  a  uniformly  illuminated 
diffuse  surface.  If  the  diameter  of  the  lens  is  d  and  the 
image  is  formed  a  distance  z  from  the  lens,  b  »  1.5  ,  as 

5  D 

shown  in  equation  (124),  while  a  more  rigorous  analysis  leads 
to  a  very  similar  expression 

bs  -  1.22($  (125) 

See  Fig.  (11)  for  the  Rayleigh  criterion  analysis  of  resolution. 

If  the  imaging  system  is  focused  on  a  relatively 
distant  plane,  z  *  f  *  focal  length  of  the  lens;  therefore 

»s  ■  i-22(4)  (126) 

where  f/D  is  the  F  number  of  the  lens;  typical  imagings 
systems  vary  from  about  f/1.4  to  f/32.  See  equation  (148) 
for  a  more  exact  equation  for  the  speckle  size. 

The  speckle  pattern  formed  by  an  imaging  scattered 
He  Ne  laser  light  (  A  *  632.8  nm)  produces  the  corresponding 
speckle  size  which  varies  from  1  to  24  micron.  See  Table  3. 
Figures  18 (a  -  f)  shows  the  correlation  between  F  number  and 
speckle  size. 

For  the  experimentation  of  this  thesis,  a  diffusely 
reflecting  object  was  viewed  through  a  lens  or  imaging  system, 
as  indicated  in  Fig.  19a. 

2.  Speckle  Photography 

As  was  stated  in  the  introduction  to  speckle  metrology. 


«***«*er  diffused  laser  light  is  passed  through  a  finite 


aperture,  a  speckle  pattern  is  formed.  For  the  science  of 
holography,  speckle  is  considered  noise  because  of  its 
limiting  effect  on  fringe  clarity.  /~  1,  217 

In  Fig.  19d,  a  double-exposure  soecklegram  is  seen 
at  a  magnification  of  30 OX.  Notice  the  speckle  pairs.  These 
speckle  pairs  developed  due  to  the  correlated  motion  of  the 
actual  speckle  pattern  during  the  double -exposure  recording. 
Processing  of  this  double -exposure  specklegram  produces  a 
Young's  fringe  pattern  from  which  in-plane  translation  and 
rotation,  strain,  and  vibrational  analysis  can  be  analyzed. 

There  are  two  prominent  advantages  of  the  speckle 
process  over  holography.  These  advantages  are  simolicitv  of 
the  optical  system  required  for  both  recording  and  processing 
the  image  and  the  relative  ease  of  interpreting  the  obtained 
results. 


- *| 

I  mag*  (film) 
plan* 

(a)  Photographic  plate  imaging  system 

Figure  19.  Two -exposure  speckle  photography  for  measure¬ 
ment  of  in-plane  translation. 
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(b)  Speckle  transparency  processing  system 


(c)  Fringe  pattern  formed  in  the  image  Diane  of  recording 
system.  The  object  was  translated  in  the  x  direction 
between  exposures. 


Figure  19.  (Continued) 
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(d)  Double- exposure  specklegram  magnified  3C0X. 

Figure  19.  (Continued) 

In-plane  displacement  will  be  measured  for  both 
translation  and  rotation  of  a  rigid,  macroscopically-f lat 
object  using  the  experimental  setup  described  by  Figure  19a. 

An  image  of  the  object  surface  is  formed  in  the  film  plane 
by  a  lens  of  focal  length  f  and  a  diameter  D.  The  object 
distance  4  and  the  image  distance  4  j  are  related  by  the 
thin  lens  imagery  equation: 

I+I  1 

l.  l>  f  (127) 

This  equation  is  approximately  valid  for  real  lenses  and  is 
quite  accurate  when  an  image  is  formed  using  only  rays  which 
are  close  to  the  optical  axis. 

Imagery  systems  or  lenses  are  specified  in  terms  of 
their  relative  aperture  or  F  number.  The  F  number  is  the  ratio 


of  the  focal  length  f  to  the  diameter  of  the  aperture  D.  If 
the  same  distant  object  is  imaged  by  a  lens  of  equal  F 
number,  the  average  irradiance  in  their  image  planes  will  be 
identical. 

For  an  object  illuminated  by  incoming  laser  light  the 
image  formed  in  the  film  plane  will  be  modulated  by  a  random 
speckle  pattern  having  a  characteristic  size  b  determined  by 
the  aperture  of  the  imaging  lens,  as  is  seen  in  equation  (126). 

3.  In-olane  Translation  Measurements  (x  and  y  Trans¬ 
lation  and  Rotation) 

If  an  object  translates  horizontally  by  an  amount 
Lx  the  relative  phase  of  the  light  in  each  of  the  various 
rays  that  contributes  to  the  formation  of  the  speckle  pattern 
will  be  unchanged.  The  entire  speckle  pattern  will  translate 
in  the  film  plane  of  the  imaging  lens  for  small  translations 
of  the  object  to  insure  there  is  no  decorrelation  of  the 
speckle  pattern.  The  soeckle  pattern  moves  a  distance  MLx  in 
the  image  plane  if  the  object  is  translated  on  amount  Lx,  where 
M  is  the  magnification  of  the  imaging  lens.  Similarly  the 
speckle  pattern  in  the  object  plane  will  translate  an  amount 
MLy  if  the  object  is  translated  an  amount  L^. 

The  amount  of  speckle  motion  is  independent  of  the 
angle  of  illumination  for  in-plane  translation. 

To  measure  in-plane  translation  a  double  exposure 
is  made.  The  film  is  exposed  once  prior  to  the  object  motion 
and  once  after  the  object  motion.  If  the  translation  is 
greater  than  the  speckle  size  b_  then  the  resultant  negative 
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of  the  image  plane  will  contain  pairs  of  identical  speckle 
patterns  as  is  seen  in  Fig.  19d.  The  speckle  pairs  will  be 
separated  by  a  distance  ML.  The  separation  ML  of  the  speckle 
pair  can  be  measured  directly  by  microscopic  examination  of 
the  film.  Alternatively,  the  translation  can  be  measured  by 
the  optical  system  seen  in  Fig.  19a.  For  this  Fourier  trans¬ 
forming  system  the  beam  is  first  collimated.  It  is  then  sent 
through  the  transforming  lens.  The  specklegram  is  then 
placed  behind  the  Fourier  transforming  lens.  The  collimated 
beam  then  passes  through  the  transforming  lens  and  the 
specklegram  and  the  speckle  pairs  recorded  on  the  specklegram 
form  Young's  fringes  in  the  back  focal  plane  of  the  trans¬ 
forming  lens  as  is  seen  if  Fig.  19c.  The  bright  central  spot 
seen  in  Fig.  19c  is  formed  by  undiffracted  light  transmitted 
by  the  specklegram.  Each  pair  of  the  corresponding  speckles 
act  as  a  pair  of  identical  sources  of  coherent  light  which 
form  Young's  fringes.  Since  all  speckle  pairs  are  separated 
by  the  same  distance  ML,  all  of  the  Young's  fringes  overlay 
to  form  the  pattern  shown  in  Fig.  19c. 

Interpretation  of  the  fringes  is  straight  forward. 

As  in  a  Fourier  transform,  the  fringe  orientation  is  normal 
to  the  in-plane  translation  L.  The  magnitude  L  of  the  trans¬ 
lation  can  be  determined  by  applying  equation  (117).  The 
speckles  on  the  transDarency  are  separated  by  a  distance  d  ; 

xfT 

the  fringes  will  have  a  soacing  df  *  jp—  .  The  in-plane 

s 

translation  of  the  object  is  given  therefore  by 
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(128) 


where  A  is  the  wave  length  of  the  laser  light  used  to  form 
the  fringes,  fT  is  the  focal  length  of  the  tranforming  lens, 

M  is  the  magnification  of  the  imaging  system  used  to  record 
the  specklegram,  and  d^  is  the  fringe  spacing. 

In-plane  rotation  is  the  clockwise  or  counter¬ 
clockwise  rotation  of  a  rigid  body.  The  recording  or 
imaging  system  for  in-plane  rotation  is  identical  to  that  for 
in-plane  translation.  Another  type  of  optical  processing 
system  must  be  utilized.  Due  to  the  nature  of  rotational  rigid 
body  motion  every  point  on  the  body  will  undergo  a  differing 
amount  of  translation.  Therefore,  the  specklegram  must  be 
analyzed  in  a  point-by-point  fashion.  The  optical  processing 
system  for  in-plane  rotation  is  seen  in  Fig.  22a.  Figure  22b 

i 

shows  the  different  fringe  patterns  that  are  obtained  from  a 
counter-closkwise  rotation  about  the  plates  center.  Th^  fringe 
spacing  is  again  inversly  proportional  to  the  point  motion  while 
the  fringes  are  90°  to  the  direction  of  the  relative  motion 
of  the  body. 

4.  Out-Of-Plane  Rotation  Measurement  (TILT)  /1T2,  53 7 

Tilt  is  out-of-plane  rotation.  Tilt  will  be  measured 
using  a  two-exposure  speckle  photograph.  The  experimental 
setup  is  identical  to  that  seen  in  Fig.  19a  utilized  for 
in-plane  translation  with  a  slight  modification.  Instead  of 
focusing  the  imaging  system  on  the  image,  the  focus  will  be 
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placed  in  the  back  focal  plane  or  at  infinity.  Again,  a 
speckle  pattern  with  characteristic  size  bg  =  1.22XF  will 
be  formed  in  this  plane. 


y 


Figure  20.  Two-exposure  speckle  photography  for  measurement 
of  tilt.  Light  rays  contributing  to  the  forma¬ 
tion  of  one  speckle  are  shown.  From  Vest  /  53 __7 

Each  speckle  is  formed  by  all  light  scattered  in  a 
Darticular  direction,  as  is  indicated  in  Fig.  20.  If  the 
object  is  translated  in  its  own  plane  L  =  L  £  +  L  j,  the 

^  y 

speckle  pattern  remains  stationary;  however  if  the  object 
tilts  by  an  angle  y  ,  as  in  Fig.  20,  the  speckle  pattern 
in  the  back  focal  plane  will  translate.  Each  ray  contributing 
to  the  formation  of  a  single  speckle  leaves  the  object  surface 
at  a  different  location  y,  but  all  travel  in  nearly  the  iden¬ 
tical  direction.  When  the  object  is  tilted  by  an  amount  y, 
there  is  a  systematic  change  of  relative  phase  among  these 
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rays.  The  phase  shift  can  be  written 


A*-^(k2-k,)*L, 


(129) 


where  K^is  the  unit  direction  vector  of  the  illumination, 
and  Kj  is  the  unit  direction  vector  of  the  rays  contributing 
to  the  speckle  under  consideration.  For  the  configuration  of 
Fig.  20,  L  *  -yyK  and 


ku  ”  0,  kiy  -  -sin ku  -  -cos0„ 
and  ku-0>  **»1; 

therefore 


(130) 

(131) 


(132) 


Equation  (132)  describes  a  linear  variation  of  phase 
with  y,  so  the  effect  is  as  if  all  the  rays  contributing  to 
a  single  speckle  tilted  by  a  small  angle,  (1  +  cosO^y  . 

The  corresponding  displacement  in  the  focal  plane  can  be  cal¬ 
culated  using  equation  (133). 


Figure  21.  Spatial  filtering  light  of  an  off  axis  plane 

wave  focused  to  an  off  axis  spot  by  a  thin  lens. 
From  Vest  /” 5^7 
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Figure  21  shows  what  can  be  considered  to  be  a  plane 
wave  indicent  from  an  off-axis  object  point  located  a  large 
distance  away,  that  is,  l  ■+  00  .  This  light  will  be  focused 

to  an  off-axis  spot  in  the  back  focal  plane.  If  the  plane 

wave  propagates  at  an  angle  with  respect  to  the  optical  axis 
it  will  focus  to  a  point  a  distance  d^  from  the  axis,  where 

df  =  f  tana  (133) 

We  are  concerned  primarily  with  the  paraxial  case; 
hence  a  will  be  small,  and 

tan  a  =»  sin  a  =  a  (134) 

Then  the  translation  of  the  speckle  in  the  back  focal  plane, 
d  ,  will  be 

5 

d  *  f(l  +  cos  9.  )Y  (135) 

S  X 

where  fT  is  the  focal  length  of  the  transformincr  lens.  If 
the  angular  aperture  of  the  lens  is  not  too  large,  the  trans 
lation  of  the  entire  speckle  pattern  will  be  given  by 
equation  (  335) . 

To  measure  tilt,  the  film  is  placed  in  the  back 
focal  plane  of  the  lens.  If  a  35  mm.  camera  is  used,  the 
camera  is  focused  at  infinity  instead  of  on  the  object.  A 
double -exposure  speckle  photograph  is  then  taken  of  the  ob¬ 
ject  in  the  infinite  focus  position.  The  object  is  exposed 
once  before  tilting  the  object  and  once  after  the  object  is 
tilted.  The  transparency  formed  by  develooing  the  film  is 
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then  illuminated  by  a  converging  spherical  wave  formed  by 
a  lens  of  focal  length  fT,  as  in  Figure  19b.  This  will  form 
a  speckle  pattern  modulated  by  cosinusoidal  fringes  of  spacing 
dj  »  i  fT/dg  *  SO  the  of  ob3ect  is  given  by 


/(!  + cos 


(136) 


Several  extensions  and  applications  of  the  pre¬ 
viously  discussed  concepts  of  speckle  photography  exist  for 
the  analysis  of  vibration  and  strain.  Theory  for  each  of 
these  two  areas  will  be  discussed  here  briefly. 

5.  Vibrational  Analysis  /~6,  20,  30,  51,  53/ 

In-plane  vibrations  can  be  measured  by  time-averaqe 
speckle  photography.  /”6,  20 7  If  the  object  of  Fig.  19a 

vibrates  harmonically  as  a  rigid  body,  L  (t)  =*  (Lx.  +  Ly . )  cos  wt 

t  1  3 

2  2  ” 

Where  L  *  (Lx^  +  Ly j  )  is  the  amplitude  of  the  vibration. 

A  time-average  exposure  will  record  speckle  streaks  of  length 
2 ML.  The  3peckles  spend  most  of  their  vibrational  period 
near  two  ends  of  the  streaks;  therefore,  if  the  film  is  pro¬ 
cessed  using  the  system  of  Figure  20b,  the  fringes  observed 
in  the  back  focal  plane  will  be  similar  to  those  caused  by  an 
object  translation  of  2L.  Specifically,  the  form  of  the 
fringes  that  modulate  the  speckle  pattern  in  the  back  focal 
plane  is  described  by  the  square  of  the  zero  order  Bessel 


Function 


where  x'  is  the  coordinate  parallel  to  the  direction  of 
motion.  The  amplitide  of  the  vibration  is  determined  by 
measuring  the  distance  df  between  the  dark  fringes  to  the 
left  and  those  to  the  right  of  the  bright  spot.  The  fringes 
occur  where  the  argument  of  JQ  equals  ±  2.4,  so 

L-0J6(%S}  <138) 

where  A  is  the  wave  length  of  the  laser  light  used  to  form 
the  fringes,  fT  is  the  focal  length  of  the  fringe  imaging 
lens  (transforming  lens),  and  M  is  the  magnification  of  the 
imaging  system  U3ed  to  record  the  speckle  photograph. 

6.  Strain  Analysis  /T,  4,  9,  15,  31-34,  40,  527 

If  the  object  under  study  is  now  strained  during  the 
experiment,  displacement  varies  from  point  to  point  on  the 
object  surface.  Hence  a  different  method  of  analyzing  the 
speckle  photographs  must  be  devised. 

Using  the  system  of  Fig.  22,  the  speckle  photograph 
can  be  examined  one  region  at  a  time  by  illuminating  it  with 
a  small  beam  such  as  an  unexpanded  He-Ne  laser  bean  as  shown 
in  Fig.  22a.  The  speckle  pattern  on  the  transparency  will 
diffract  this  beam  into  a  diverging  cone  of  angular  extent 


where  f/D  is  the  F  number  of  the  lens  used  to  record  the 
speckle  photograph. 


(b)  Appearance  of  fringes  at  various  locations  on  the 
specklegram  of  a  rigid  plate  that  was  rotated  about 
a  normal  axis  through  the  center  0. 

Figure  22.  Formation  of  fringes  by  illuminating  a  double 
exposure  photograph  with  a  thin  laser  beam. 

From  Vest  /“ 53 __7 

If  the  small  illuminated  region  contains  oairs  of 

identical  speckles  displaced  by  a  distance  dg,  the  light  in 

the  resulting  pairs  of  diffraction  cones  will  interfere  to 

form  a  diffraction  pattern  modulated  by  Young’s  fringes.  ^7#  4, 

537  If  the  fringes  are  viewed  from  the  photograph  at  a  dis- 

Xz 

tance  Z  from  the  transparency#  their  spacing  will  be  df  =  g —  . 

s 

Since  dg  ■  ML#  where  M  is  the  lateral  magnification  of  the 
optical  imaging  system  used  to  record  the  speckle  Dhotograph 
and  L  is  the  unknown  displacement# 


L 
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(140) 


Equation  (140)  gives  the  displacement  of  points  in 


a  single  small  region  of  the  object  surface.  The  fringes 
are  normal  to  the  local  in-plane  motion. 

An  alternative  method  of  analyzing  strain  was  dis¬ 
covered  by  Professor  Y.  Y.  Hung.  It  is  the  Speckle  Shearing 
Camera  /” 32 J.  This  method  utilizes  an  image  shearing  inter¬ 
ferometer  camera. 


COHERENT 
UGH  I  SOURCE 


IMAGE  SHEARING 
CAMERA 

(a)  Schematic  diagram  of  optical  setup  of  an  image 
shearing  camera. 


HEDGE 


(b)  Imaging  detail  of  the  image  shearing  camera. 


Figure  23.  Two  exposure  speckle  photography  using  image 
shearing  camera  for  strain  and  vibrational 
measurements.  From  Hung  /“32 J 

The  image  shearing  camera  of  Fig.  23a  has  its  lens 

divided  into  two  parts  by  a  thin  glass  wedge.  The  lens 

images  the  object,  which  is  illuminated  by  a  diverging  beam, 
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inclined  at  an  angle  0^  to  the  viewing  direction.  Since  the 
rays  passing  through  the  glass  wedge  are  deviated,  the  two 
images,  (one  focused  by  each  half  the  lens)  are  laterally 
sheared  with  respect  to  each  other  as  is  seen  in  Fig.  23b* 

Let  the  plane  containing  the  prism  angle  be  parallel 
to  the  x-z  plane.  Rays  passing  through  the  wedge  are  laterally 
shifted  in  the  x-direction.  Assuming  that  the  wedge  is  placed 
very  close  to  the  lens,  the  amount  of  shift  dx'  in  the  image 
plane  is 

ix'  ■  H,U i  —  1)<«. 

<141> 

Where  is  the  image  distance  from  the  lens,  M  is  the  re¬ 
fractive  index  of  the  prism,  and  a  is  the  prism  angle.  The 
equivalent  shift  5x  on  the  object  is 


4*  ■ 


(142) 


Where  M  is  the  magnification,  which  is  the  ratio  of  the 
image  distance  to  the  object  distance  Dq  from  the  lens. 
Thus, 


5x  *  fix’  <Do/Di*  *  Do*y  ”  11  a 

The  image  shearing  camera  brings  the  rays  from  point 
P(x,y)  on  the  object  to  meet  with  that  from  the  neighboring 
point  P(x  ♦  dx,y)  in  the  image  Diane.  Since  the  object  is 
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coherently  illuminated,  the  rays  from  the  two  points  inter¬ 
fere  with  each  other  in  the  film  plane. 

The  image  shearing  camera  focuses  two  laterally- 
sheared  coherent  images  in  the  film  plane  which  interfere 
with  each  other.  The  direction  of  shearing  depends  on  the 
orientation  of  the  plane  containing  the  wedge  angle.  Thus, 
by  rotating  the  wedge  on  the  wedge-lens  assembly  about  the 
axis  parallel  to  the  viewing  direction,  the  direction  of 
shearing  can  be  varied.  The  amount  of  shearing  can  be 
adjusted,  as  indicated  in  equation  (143) »  by  varying  the  wedge 
angle  a  as  well  as  the  object  distance  DQ.  When  the  object 
is  deformed,  a  relative  displacement  between  two  points  pro¬ 
duces  a  relative  optical  phase  change  A  given  by: 


^  j(l  +  COS0  )  [*(x4dx,y)  -w(x,y)j 
sine[u(x+dx,y)  -  a(x,y)]J 


(144) 


where  u  and  w  are  the  displacement  components  in  the  x  and 
z  directions,  respectively*  If  the  shear  dx  is  small,  the 
relative  displacemnts  may  be  approximated  by  the  disolace- 
ment  derivatives  and  thus  equation  (144)  becomes 


A 


(145) 


where  the  strains  are  given  by 


c 


x 


3u 

lx 


3v 


3v 

Tx 
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The  term  fix  is  equivalent  to  the  gage  length.  Should  the 
shear  be  in  the  y  direction,  the  derivatives  in  equation  (145) 
would  be  with  respect  to  y.  This  is  achieved  by  rotating  the 
camera  lens  90°.  If  the  object  is  rotated  90°  about  the  z 
axis,  u  is  replaced  by  v,  the  displacement  component  in  the 
y  direction. 

By  taking  a  double-exposure  recording  with  the  object 
being  strained  between  exposures,  a  speckle  fringe  pattern 
depicting  A  will  be  generated.  Dark  fringes  occur  when 

A  -  Ntt  (146) 

where  N  »  1,  3,  5,  7  ... 

3w 

To  measure  this  derivative  may  be  isolated 
using  normal  illumination  (6  »  0). 

However,  it  is  impossible  to  isolate  3u/3x.  To 

3u 

achieve  this  isolation  requires  recording  two  fringe 
patterns  using  two  different  illumination  angles;  then 

3w 

can  be  separated  from  -g-j  on  a  point  by  point  basis  using 
the  following  equation: 


du 

Tx 


X 

77x 


N^(l  +  cos62)  -  N2(l  +  cos0^) 
sinO^  (1  +  cos02)  -  sin02  (1  +  cos0^). 


(147) 


where  and  N2  are  fringe  orders  corresponding  to  illumi¬ 
nation  angles  0^  and  02  respectively. 
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G.  FILM  THEORY 


1.  Film  And  Optical  Requirements 

The  sensitivity  measurement  by  speckle  Dhotography 
is  limited  by  the  requirement  that  speckle  motion  in  the  film 
plane  must  exceed  the  characteristic  speckle  size, 

if 

b  ■  1.22  jj—  .  If  f/2.8  is  used  as  a  camera  aperture, 
then  bg  ■  2.16pm  and  a  film  resolution  of  about  462  lines/mm 
is  required.  A  high  resolution  film  is  necessary  for  speckle 
photography.  S0253  film  was  used  for  all  SDeckle  experimen¬ 
tation.  This  film  has  a  resolution  of  2500  lines/mm.  The 

_ J _ _ _ _ _ - _  _  c  t _ _ l  _ _ i.  j _  ■  _  e  /  f\  _  r 


maximum  frequency  of  the  speckle  pattern  is  f. 


B) 


564. 


Since  film  speed  varies  inversely  with  resolution,  relatively 
long  exposure  times  were  required. 

2.  Measurement  Limitations 


The  upper  limit  to  measurable  displacements  or  tilts 
is  set  by  the  decorrelation  of  the  two  speckle  patterns.  This 
is  evident  because  processing  of  speckle  photographs  requires 
formation  of  Young's  fringes  by  diffraction  of  light  by  nairs 
of  identical  speckles.  A  loss  of  correlation  will  be  caused 
if  strains  are  large  enough  to  alter  the  microscopic  structure 
of  the  opaque  surface  under  study.  Loss  of  correlation  also 
occurs  when  motions  are  sufficiently  large  that  the  aperture 
of  the  object  forming  lens  samoles  an  appreciably  different 
portion  of  the  wave  front  scattered  by  the  object  during 
the  two-exposures.  When  this  occurs,  only  a  fraction  of  the 
speckle  pattern  recorded  on  the  photograph  consists  of  oairs 
of  identical  speckles,  so  the  fringes  they  generate  are  washed 
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out  by  the  light  diffracted  by  the  rest  of  the  speckle.  A 
rough  estimate  of  the  sensitivity  range  of  speckle  photography, 
based  on  experiments  reported,  is  that  motions  resulting  in 
speckle  translations  of  the  order  of  15-200  ym  in  the  film 
plane  can  be  measured.  /~ 6 ,  20,  5^7  In  the  upper  limit,  the 
object  translation  equals  about  ten  percent  of  the  diameter 
of  the  entrance  pupil  of  a  typical  camera  lens  at  f/4. 

The  practical  sensitivity  range  in  terms  of  actual 
object  translation  or  tilt  depends  on  object  size,  magnification, 
and  required  depth  of  focus.  £""21,  447  In  general,  the  range 
extends  to  motions  as  large  as  about  1.5  mm.  Above  this  limit, 
direct  measure  by  ordinary  double -exposure  photography  is  more 
practical.  The  lower  limit  of  the  range  is  about  5ym  ,  but 
depth  of  focus  requirements  usually  raise  this  limit  to  about 
100  yrn  or  higher.  Hence,  in  practice  there  is  a  large  gap 
between  the  maximum  displacement  measurable  by  holographic 
interferometry  and  the  smallest  displacement  measurable  by 
speckle  photography.  /—44  7 
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III.  EXPERIMENTAL  ANALYSIS  AND  PROCEDURES 


Two  types  of  film  were  used  during  this  thesis  experi¬ 
mentation.  During  the  Fourier  transform,  speckle  micro¬ 
photography,  and  basic  Fraunhofer  diffraction  studies,  Polaroid 
positive  negative  4x5  land  film  with  an  ASA  (55)  was  used. 

An  f/16  aperture  and  an  exposure  time  of  between  1/250  sec. 
and  1/125  sec.  were  used  for  optimim  results.  A  Kodak  special 
order  (S0253)  film  was  used  for  all  speckle  photography.  A 
Pentax  35  mm  camera  was  used  with  either  a  55  or  85  mm  focal 
length  lens.  For  object  speckle  photography  using  the  85  mm 
lens  with  an  M  *  .11  and  an  F  ■  2.8  for  optimum  results  an 
exposure  time  of  .5  sec.  was  used.  For  object  speckle  photo¬ 
graphy  using  the  55  mm  lens  with  an  M  ■  .137  and  an  F  *  4 
and  an  object  distance  of  450  mm  an  exposure  time  of  2  seconds 
was  used  for  optimum  results.  For  fringe  photography  only 
the  55  mm  lens  was  used  with  an  F  ■  2.0  with  an  object  distance 
of  five  feet  giving  exposure  times  of  5-7  min.  depending  on 
the  fringe  spacing. 

The  resolution  of  the  S0253  film  is  very  high  (on  the 
order  of  2500  lines/mm) .  The  higher  resolution  demanded 
longer  exposure  times,  but  this  higher  resolution  was  needed 
as  will  be  seen  during  the  discussion  of  speckle  size. 

Lastly,  the  film  developing  was  done  in  Kodak  D-19  de¬ 
veloper  70°F  for  5*j  minutes  with  constant  agitation.  Then 
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using  Kodak  rapid  fixer/  the  film  was  fixed  for  3  minutes. 

A  film  bath  followed  which  usually  required  30-40  minutes 
of  constant  water  flow  due  to  the  heavy  blue  tint  of  the 
S02S3  film. 

B.  SPECKLE  SIZE 

As  was  seen  in  the  theory  section  under  speckle  size 
from  the  governing  equation  (126)  speckle  size  is  a  function 
of  the  wave  length  of  light  and  the  F  number  of  the  imaging 
lens.  Table  3  tests  the  various  F  numbers  and  their  corre¬ 
sponding  speckle  size/  b  ,  when  using  a  He-Ne  gas  laser  as 

5 

was  used  in  the  experimentation  for  this  thesis.  The  wave 
length  of  a  He-Ne  gas  laser  is  A  »  . 6328. 10“3mm. 


TABLE  3 


Speckle  Size  Versus  Lens  F  Number  And  Film  Resolution 


F  1.9 

2  2.8  4  5.6  8 

11 

16 

32 

D  „x  (10~3mn$  1.4668 

1.544  2.16  3.088  4.323  6.176 

8.492 

12.352 

24.704 

resolution 

(lines/mm) 

647  462  323  231  161 

117 

80 

Looking  at  the  optical  system  used  to  record  the  speckle  , 
lens  magnification  should  be  included  in  the  speckle  calcu¬ 
lation.  /~3/  47  Fig.  -24  shows  the  ob ject-to-image  oosition 
which  relates  the  lens  equation  to  the  magnification. 
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Figure  24.  Formation  of  a  real  image  with  a  thin,  spherical 
lens. 

Assuming  the  size  of  the  speckle  dot  to  be  equal  to  the 
Airy  disk  in  the  image  plane 

bg  *  from  equation  (125). 

where  a  *  z  in  thi3  case.  From  the  lens  equation  (127) 


f  =*  or  f(a  +  b)  *  ab 

which  relates  the  speckle  size  to  the  image  magnification, 
M.  It  is  seen  that 

K  -  1.22Xf  (a  +  b)  _  1. 22Xf  /,  .  a \ 


MAGNIFICATION  FACTOR 


v/here 

M  *»  — g  is  the  magnification  of  the  imaging  system. 

Fig.  25  shows  the  relationship  between  the  magnification 
M  and  the  speckle  size  b  . 

S 

Another  area  of  concern  is  the  effect  of  magnification 
change  during  laser  speckle  photography.  /”45 7  An  assumption 
is  made  that  during  the  speckle  photograph  no  magnification 
change  occurs#  where  in  fact  a  change  does  take  place. 

For  example,  in  the  plane  stress  case,  a  magnification 

change  arises  from  the  Poisson's  effect  whereby  the  plane 
model  undergoes  lateral  thickness  changes  when  loaded  in 

its  plane. 

In  the  tilt  case,  a  magnification  change  occurs  due  to 
the  out-of-plane  motion  toward  the  camera  thus  changing  the 
object  distance  and  the  magnification  M.  This  effect  of 
change  of  magnification  is  examined  in  detail  by  Archbold 
and  Ennos  /~ 5 ,  7,  1 £7  and  Dedretti  and  Chiang  /— 457  The 
problem  of  magnification  change  can  be  minimized  by  usinq 
imaging  lenses  of  long  focal  lengths. 

The  Airy  disk  size  limitation  on  speckle  size  consider¬ 
ations  for  resolution  is  a  good  approximation,  but  subpeckle 
size  measurements  can  be  effected.  /~*54_7  This  lower  measure¬ 
ment  range  can  easily  be  reduced  to  one  order  of  magnitude  of 
the  usual  Airy  disk  limit.  The  method  of  subspeckle  size 
speckle  photography  measurement  is  discussed  in  detail  by 
Vikram  and  Vedam  £3’4_7.  It  requires  a  known  lateral  movement 
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being  measured.  The  orientation  of  the  Young's  fringes 
gives  the  unknown  deformation. 

There  are  two  other  fundamental  limitations  in  the 
recording  methods  of  this  thesis.  First,  laser  speckles  are 
three  dimensional.  Speckle  patterns  decorrelate  in  a  distance 
which  is  small  compared  to  conventional  depth  of  field  for 
an  imaging  lens  system  /~3,  4 7.  Decorrelation  is  particularly 
severe  with  lens  apertures  wider  than  about  f/4.  Most 
speckle  photography  of  this  project  was  done  with  a  lens 
aperture  of  f/2.8  therefore  decorrelation  of  speckle  was  a 
problem.  An  increased  image  magnification  would  reduce  the 
problem  of  the  lens  aperture  effect,  but  again  a  low  magni¬ 
fication  was  used  for  most  experimentation  (M  *  .11). 

Second,  speckle  photography  is  significantly  vulnerable  to 
lens  aberrations.  Ennos  recognized  this  effect  and  published 
some  data  with  regard  to  in  plane  rotation.  /~5,  107.  A  more 
indepth  study  of  the  problem  is  given  by  Stetson  /“’46_7. 
Stetson's  figures  showing  a  comparison  of  Young's  fringes 
achieved  by  a  well-corrected  copy  lens  at  unity  magnification 
and  speckle  halos  from  a  normal  camera  at  the  same  magnifica¬ 
tion  vividly  illustrate  the  loss  of  fringe  visibility  due  to 
lens  aberations. 

C.  FOURIER  TRANSFORMS 

In  order  to  better  understand  diffraction,  lens  resolution 
spatial  filtering  and  other  areas  of  optical  physics,  a  pre¬ 
liminary  study  of  Fourier  transforms  was  conducted. 


The  first  experiment  that  was  conducted  consisted  of 
comparing  differing  rectangular  aperture  spacing.  Fourier 
transforms  were  achieved  using  Ronchi  ruling  combinations  of 
100-300  lines/mm  and  100-450  lines/mm.  The  results  of  the 
100-300  lines/mm  Ronchi  ruling  is  seen  in  Fig.  26. 


Figure  26.  Fourier  transform  of  a  100  and  300  lines/mm 
Ronchi  ruling  combination. 

Three  types  of  filtering  were  conducted.  Using  two 
Ronchi  rulings  of  300  lines/mm,  a  rectangular  aperture  grid 
was  setup.  First,  the  center  of  the  Fourier  transform  of  the 
rectangular  aperture  was  filtered,  which  eliminated  all  lines 
from  the  reproduced  grid  pattern  as  seen  in  Fig.  27a. 

Second  ,  the  horizontal  lines  were  filtered  out  of  the 
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Fourier  transform  of  the  rectangular  aperture,  which  elimi¬ 
nated  all  horizontal  lines  from  the  reproduced  grid  pattern. 
Last  ,  the  vertical  lines  were  filtered  out  of  the  Fourier 
transform  of  the  rectangulary  aperture,  which  eliminated  the 
vertical  lines  from  the  grid  pattern  as  is  shown  in  Fig.  27b. 


(a)  Image  pattern  as  a  result  of  filtering  the  center 
of  the  Fourier  transform. 


(b)  Image  pattern  as  a  result  of  filtering  the 
vertical  lines  of  the  Fourier  transform. 

Figure  27.  Fourier  transform  filtering  using  a  rectangular 
grid. 

Figure  28  shows  the  lab  setup  for  all  Fourier  transform 
experimental  analysis.  A  1.5  milliwatt  He-Ne  gas  laser  with 
a  *  *  632,8  nm  was  used  to  produce  the  monochromatic  wave. 

A  combination  10X  microscope  objective  combined  with  a  25  iim 
spatial  filter  were  used  to  expand  the  beam.  All  lenses  in 
this  Fourier  transform  arrangement  were  spherical  double  con¬ 
vex  lenses  with  a  focal  length  fT  of  177.8  mm.  Another  lens 
collimated  the  beam.  The  grid  was  then  placed  in  the  optical 
system  followed  by  the  transforming  lens.  The  grid  pattern 
was  then  reproduced  after  passing  through  the  beam  filter. 
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Figure  28.  Equipment  arrangement  for  conducting  Fourier  transform  experimentation. 
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A  10  X  microscone  lens  was  used  to  expand  the  beam  onto  a 
mirror  which  then  reflected  the  reproduced  grid  image  on  a 
frosted  glass  screen.  A  magnifying  eyepiece  enabled  the 
observer  to  view  the  imaqe  grid. 

In  order  to  process  the  double -exposure  specklegram,  the 
Fourier  transform  system  of  Fig.  28  was  utilized.  This  system 
of  fringe  analysis  was  done  using  either  the  film  prior  to 
or  after  the  transforming  lens  as  is  seen  in  Fig.  15  b  &  c 
and  described  in  detail  in  the  Fourier  transform  section  of 
this  thesis. 

D.  IN-PLANE  TRANSLATION  AND  ROTATION 

Figure  29  shows  the  lab  setup  for  in-plane  translation 
and  rotation.  A  1.5  milliwatt  Ke-Ne  gas  laser  was  used  to 
produce  a  monochromatic  wave. 

A  double  shutter  was  used  to  facilitate  the  double -ex¬ 
posure  speckle  Dhotographing  system.  A  40x  lens  with  a  25  um 
spatial  filter  was  used  to  exoand  the  beam.  A  5"  x  7" 

2024-T4  Aluminum  plate  was  used  for  the  object.  A  Pentax 
35  mm  camera  with  an  85  mm  lens  was  used  to  record  the  image 
at  an  object  distance  of  850  mm.  Optimum  exposure  times  for 
this  setup  were  second  for  an  F  ■  2.8. 

For  in  plane  translation,  a  micrometric  holder  which  read 
in  1/10,000  in.  and  1  turn  ■  1/40  in.  was  used  while  for 
in-plane  rotation  a  micrometer  holder  which  read  in  .01  degrees 
and  1  turn  >0.1  degrees  was  used. 

Translation  of  .003  in.  to  .070  'in.  were  measured. 

A  plot  of  actual  translation  versus  fringe  calculated 
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Figure  29.  Laboratory  arrangement  for  recording  in-olane  translation  and  rotation 


translation  can  be  seen  in  Fig.  30,  while  Table  4  shows 
calculation  of  fringe  spacing  and  fringe  density  for  each 
translation. 

The  average  error  for  the  overall  translation  data  was 
calculated  to  be  one  percent.  Both  Fig.  30  and  the  calcu¬ 
lated  error  show  that  there  was  excellent  correlation  between 
the  actual  and  calculated  translation.  Fig.  31  (a-f)  show 
fringe  spacing  for  various  translations.  As  in  the  Fourier 
transform  experimentation  done  previously,  it  can  be  seen 
that  for  small  translations,  speckle  pair  spacing  is  small, 
the  aperture  of  diffraction  on  the  specklegram  is  small  and 
the  resultant  Young's  fringe  spacings  are  large,  and  con¬ 
versely  for  large  translation  motion,  the  speckle  pair  spacing 
is  large  on  the  specklegram,  therefore  a  large  diffraction 
aperture  resulting  in  Young's  fringes  which  are  closely  spaced 
together. 

Figure  33  shows  fringes  produced  from  specklegrams  of 
various  rotations.  Seeing  that  each  point  in  the  specklegram 
field  is  rotated  a  different  amount,  individual  points  on  the 
specklegram  were  inspected. 

As  seen  in  Fig.  32-33,  fringes  showing  almost  pure  trans¬ 
lation  in  both  the  x  and  y  direction  were  observed  at  the 
four  quarters  of  the  specklegram.  At  the  center  of  the 
specklegram  no  relative  displacement  occurs,  therefore  Young's 
fringes  are  not  visible.  At  the  outer  edge  between  the 
quarter  points  of  the  specklegram,  diagonal  fringes  show  the 
existence  of  both  x  and  y  components  of  the  rotation.  Moving 
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RCTURL  TRRNSLRTION  ;  (Mini-inches) 
Figure  30. 


TABLE  4 


1  - 

Calculation  Of 

Fringe  Spacing  For  Various  Translations 

l 

A  actual 
X0.001 
(in.  ) 

lines 

per  area 

area 

(mm) 

df 

(mm) 

L  (mm) 
exp. 

L  (in; 
X0.001 
exp 

i„.Error  *  TM 

1 

8.61 

5.5 

600 

109.09 

.2259 

8.89 

.0325 

'i 

1 

12.31 

8 

600 

75.0 

.3286 

12.9 

.0479 

- 

18.45 

12 

600 

50.0 

.493 

19.4 

.0515 

r  ■ 
h 

24.37 

16 

600 

37.5 

.656 

25.8 

.0374 

> 

30.06 

18 

600 

33.33 

.7395 

29.1 

.0319 

1 

31.98 

20 

600 

30 

.8215 

32.3 

.010 

f 

38.58 

24 

600 

25 

.986 

38.81 

.0059 

I 

40.37 
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Figure  31  (Continued) 
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(f)  Fringes  generated  due  to  0.050  inch  translation. 


Figure  31  (Continued) 
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Fiqure  32.  Double  exposure  rotation  specklegram  shov/incr 

location  of  various  fringe  observation  points. 


(a)  Lack  of  fringes  observed  at  Point  A. 


Figure  33.  Young’s  fringes  produced  from  various  points  of 
observation  on  a  rotation  specklegram  as  indi¬ 
cated  in  Figure  32. 
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(b)  Rotation  fringes  observed  at  Point  B  on 
the  specklegram. 


(d)  Rotation  fringes  observed  at  Point  D  on 
the  specklegram. 


(e)  Rotation  fringes  observed  at  Point  E  on 
the  specklegram. 


Figure  33  (Continued) 
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(f)  Rotation  fringes  observed  at  Point  F 
the  soecklegram. 


(q)  Rotation  fringes  observed  at  Point  G 
the  specklegram. 


Figure  33  (Continued) 


inward  from  each  of  these  oomts  toward  the  center  of  the 
specklegram,  the  relative  displacement  is  less,  as  would 
be  expected,  the  Young's  fringe  spacing  increases. 

E.  OUT-OF- PLANE  ROTATION  -  (TILT ) 

Figure  34  shows  the  laboratory  arrangement  of  equipment 
for  imaging  tilt  double-exposure  specklegrams. 

The  laboratory  setUD  is  identical  to  the  one  used  for 
recording  translation  and  rotation  with  two  differences. 

First,  after  the  beam  is  expanded,  it  is  collimated  prior 
to  illuminating  the  object.  This  collimation  process  is 
necessary  for  recording  all  out-of-plane  motion  using  this 
type  recording  system.  Second,  instead  of  recording  the 
image  in  the  image  plane,  the  camera  is  focused  at  infinity, 
placing  the  image  in  the  back  focal  plane  of  the  camera.  An 
F  number  of  2  was  used  at  an  object  distance  of  450  mm  with 
am  exposure  time  of  3  seconds  per  exposure  for  optimim  re¬ 
sults.  A  tilt  range  of  1.5"  -  15"  was  tested  producing  fringe 
as  seen  in  Figures  35  a  and  b. 

Calculations  were  made  of  fringe  spacing,  tilt,  and  tilt 
angle  and  compared  with  the  actual  tilt  observed  during 
experimentation.  Results  of  these  calculations  are  tabulated 
in  Table  5. 

A  graphical  comparison  of  calculated  tilt  angle  versus 
calculated  fringe  is  observed  in  Fig.  36. 

Both  the  graph,  the  tabular  results,  and  the  mean  error 
analysis  confirm  that  there  is  excellent  correlation  between 


Figure  34.  Laboratory  equipment  arrangement  for  imaging  tilt  double -exposure 
specklegrams. 
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Figure  35.  Young's  fringes  produced  from  various  angles 
of  out  of  Plane  rotation  (tilt) . 
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the  actual  tilt  ancle  observed  during  experimentation  and 
the  fringe  calculated  tilt  angle. 

P.  MEASUREMENT  OF  DISPLACEMENTS  DUE  TO  VIBRATION 

Figure  37  shows  the  laboratory  setup  for  the  in-plane 
vibrational  experimentation  that  was  conducted  in  conjuction 
with  this  thesis.  The  equipment  required  for  vibrational 
analysis  was  the  same  as  for  tilt  experimentation  with  the 
exception  of  the  plate  holders,  excitoi,  power  transformer 
and  oscilliscope. 

The  plate  holder  was  modified  to  enable  small  amplitude 
in-plane  vibrations  to  be  transmitted  to  the  5x7  aluminum 
Dlate.  A  piezoelectric  excitor  linked  in  series  with  a 
spring  having  a  large  spring  constant  enabled  vibrational 
amplitudes  of  .1  -  2  mm  range  to  be  transmitted. 

The  image  was  recorded  with  a  35mm  Pentax  camera  using 
the  double -exposure  speckle  technique  described  previously. 
An  F  ■  2.8  with  an  85mro  lens  using  an  exposure  time  of  2 
seconds  per  exposure  were  used  to  record  the  image  at  an 
image  distance  of  850  mm  and  a  magnification  of  .11. 

G.  STRAIN  ANALYSIS 

Figure  38  shows  the  experimental  setup  for  the  strain 
analysis  conducted  in  this  thesis.  With  the  exception  of 
the  plate  holder  and  a  modification  to  the  imaging  camera, 
which  will  be  discussed  later,  the  equipment  setup  is  iden¬ 
tical  to  the  tilt  experimental  equipment  system  discussed 
previously. 
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Figure  38.  Laboratory  equipment  arrangement  for  strain  analysis. 


The  image  shearing  camera  was  used  to  record  the  double 
exposure  specklegram.  In  the  theory  section,  the  image 
shearing  camera  was  described  in  general.  Both  a  55mm  and 
85mm  lens  were  used  with  the  Pentax  35mm  camera .  The  shearing 
prism  was  a  1°  glass  wedge  with  an  index  of  refraction  of  1.57. 
With  the  55mm  lens  an  image  distance  of  450mm  and  an  F  of 
2.0  with  an  exposure  time  of  3  seconds  per  exposure  was  used, 
while  with  the  85mm  lens  an  image  distance  of  850mm  and 
an  F  of  2.8  with  an  exposure  time  of  2  seconds  per  exposure 
was  used.  The  image  shearing  camera  gives  a  double-exposure 
specklegram  which  is  processed  using  the  system  seen  in  Fig.  39. 

i^iqure  40  shows  the  fringe  pattern  obtained  from  the 
processed  image  shearing  specklegram.  This  specklegram  can 
be  processed  in  white  light  without  a  lens  if  the  shearing  is 
adequate . 

To  produce  the  fringes  of  Fig.  40,  a  6  in.  x  8  in. 
rectangular  aluminum  plate  with  its  surface  sprayed  flat  red 
was  used.  The  plate  was  securely  clamped  along  its  four 
boundaries.  With  collimated  illumination  with  a  He-Ne  laser 
a  double  exposure  of  the  image  was  photographed.  The  plate 
was  displaced  at  the  center  0.0005  between  exposures. 

A  no  fringe  area  was  seen  during  experimentation  in  areas 
of  plastic  deformation.  During  plastic  deformation  there  is 
no  motion  of  that  area  of  the  body,  therefore  there  are  no 
local  fringes  developed  in  any  area  where  plastic  deformation 
is  experienced. 
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Figure  40.  Fringe  pattern  obtained  from  the  double -exposure 
specklegram  for  a  displacement  at  the  center  of 
the  plate  of  0.0005  inches. 
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Compared  with  other  methods  of  speckle  photography  and 
interferometry,  the  image  shearing  camera  was  found  to  pro¬ 


duce  superior  fringes.  Besides  the  improved  quality  of 
fringes  there  are  several  other  areas  where  the  speckle  image 
shearing  camera  improves  over  conventional  speckle  methods  or 
holography. 

First,  the  setup  is  relatively  simple  and  does  not  require 
precision  alignment  of  optics.  Second,  it  does  not  require 
the  special  vibrational  isolation  that  accompanies  holographic 
measurements.  Third,  the  need  of  coherence  in  the  illuminating 
light  source  is  minimized.  Although  a  collimated  beam  was 
utilized,  it  was  not  required.  Fourth,  it  provides  a  wider 
and  more  controllable  range  of  sensitivity.  The  sensitivity 
can  be  changed  by  either  changing  the  illumination  angle  ei 
or  the  amount  of  shearing.  Fifth,  the  strain  is  measured 
directly  by  the  image  shearing  camera. 

While  this  technique  works  well  for  small  shearing  effects, 
like  other  speckle  techniques  the  image  shearing  method  is 
highly  susceptible  to  decorrelation  when  the  object  under¬ 
goes  large,  tilt,  translation,  or  strains.  This  decorrelation 
deteriorates  the  fringe  visibility. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  following  conclusions  have  been  drawn  from  the  research 
conducted  during  the  course  of  this  project. 

A  comparison  of  holographic  interferometry  to  speckle 
metrology  shows  that  both  systems  of  measurement  have  favor¬ 
able  points  which  support  further  research  in  each  area. 

The  simplicity  of  equipment  arrangement  for  speckle 
metrology  requires  no  precision  alignment  of  optics  as  is  re¬ 
quired  in  holographic  measurements  /”7 __7*  Speckle  metrology 
requires  no  vibrational  isolation  system.  The  use  of  multi¬ 
node  lasers  is  possible  because  speckle  metrology  minimizes 
the  coherent  length  of  the  light  source  required  for  accurate 
measurements.  There  is  a  wider  range  of  sensitivity  and  an 
ease  of  sensitivity  variation  in  speckle  metrology  that  is 
not  found  in  holography  /~34 ,  35,  537.  Speckle  metrology 
requires  a  much  lower  film  resolution.  Decorrelation  of  the 
speckle  pattern  due  to  large  tilt,  translation,  or  strains 
deteriorates  the  fringe  visibility.  Speckle  produced  fringes 
are  highly  susceptible  to  lens  aberration.  Even  at  best, 
speckle  produced  fringes  are  of  poorer  quality  than  holographic 
fringes.  Holographic  interferometry  presents  a  much  broader 
band  of  measurement  /”53 _7«  Three  dimensional  motion  can  be 
studied  using  holography  much  more  readily  than  can  be 
achieved  by  speckle  metrology  7“ 53  7. 

From  the  previously  listed  conclusions,  it  is  seen  that 
holographic  interferometry  used  in  conjunction  with  speckle 
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metrology  can  provide  accurate,  wide  ranee,  highly  sensitive, 
means  of  metrology  to  Naval  Engineering  in  general  and  Naval 
shiDyards  specifically. 

Specific  conclusions  derived  by  actual  experimentation 
are  as  follows: 

Measurements  of  in-plane  translations  were  measured  for 
a  range  of  0.001  inches  to  0.07  inches  with  a  mean  error  of 
less  than  four  percent. 

Actual  measurements  of  in-plane  rotation  were  not  made, 
but  rather  a  comparison  of  relative  fringe  spacing  from  indi¬ 
vidual  point  inspection  of  the  specklegram  was  performed. 
Direction  of  rotation  and  relative  size  of  point  displacement 
could  be  determined  from  the  resultant  fringe  patterns. 

Measurements  of  out-of-plane  rotation  (tilt)  were  made. 
The  range  of  experimentation  of  1.5"  to  15"  was  too  small  and 
too  few  points  were  sampled  to  effectively  conclude  that  the 
accuracy  of  measurement  was  as  precise  as  experimentation  in¬ 
dicates.  Also  a  more  effective  method  of  measuring  actual 
tilt  should  be  used  to  insure  the  actual  tilt  measurements 
are  exact. 

Vibrational  analysis  of  both  in-plane  amplitude  measure¬ 
ments  and  modal  shape  studies  were  conducted  with  limited 
success.  The  fringes  produced  from  in-plane  vibrations  gave 
good  results,  but  too  small  an  amplitude  range  and  too  few 
points  were  sampled  to  make  any  conclusions  concerning  the 
effectiveness  of  this  system. 
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Strain  analysis  using  the  image  shearing  camera  produced 
fringes  which  showed  the  strain  field  over  the  entire  plate 
due  to  a  center  displacement  of  approximately  0.0005  inches. 
Object  and  image  distances  can  be  varied  to  increase  the 
range  of  shear  observation;  but  this  was  not  done.  Therefore, 
the  effectiveness  of  the  image  shearing  camera  for  a  wide 
range  of  measurements  can  not  be  stated.  Good  fringes  were 
achieved  using  this  system  for  a  limited  range  of  shear. 

Each  area  of  study,  i.e.,  translation  both  in-plane  and 
out-of-plane,  vibrations,  and  strain,  is  a  broad  enough  field 
sufficient  to  support  thesis  investigation. 

B.  RECOMMENDATIONS 

The  field  of  speckle  metrology  is  in  its  embryo  stage  of 
development.  New  techniques  and  methods  of  analysis  and 
measurement  are  being  discovered  continually.  The  only  re¬ 
striction  placed  on  this  field  is  the  imagination  of  the 

i 

experimenter . 

It  is  recommended  that  several  of  the  other  methods  of 
speckle  metrology  be  investigated.  Experimentation  in  two 
and  four  aperture  Fourier  filtering  analysis  of  two  and 
possibly  three  dimensions  is  being  conducted  with  considerable 
success  l]i,  7,  9,  21,  53^.  Dual  beam  speckle  interferometry 
for  strain  component  isolation  and  vibrational  analyses  is 
achieving  results  comparable  to  that  of  holography  /“15,  30, 
37,  40,  51,  537. 
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In  the  area  of  in-plane  rotation,  it  is  recommended  that 
x  and  y  component  isolation  be  conducted  to  measure  both 
direction  and  relative  amount  of  translation  using  fringe 
spacing  and  orientation. 

The  system  of  experimentation  used  to  measure  tilt  was 
found  to  be  highly  effective  and  accurate.  Two  improvements 
could  be  made  in  this  area.  First,  a  more  accurate  method 
is  needed  for  measuring  the  actual  out-of-plane  the  tilt 
which  has  occurred.  This  would  validate  the  error  analysis 
conducted  in  this  project.  Second,  an  increased  range  of 
measurement  could  be  achieved  by  varying  the  object  and 
image  distances . 

The  surface  of  vibrational  and  strain  analyses  were  barely 
scratched  during  the  course  of  this  project.  In  conjunction 
with  these  areas  of  study,  it  is  recommended  that  speckle 
metrology  be  utilized  for  beam  analysis.  Deflection,  slope, 
and  moment  measurements  could  easily  be  effected  using  the 
dual  beam  /—5, 6, 10,217,  Fourier  filtering  /~21,53 7,  or  shear 
camera  /“31-337  methods  discussed  previously.  Vibrational 
analysis  should  be  conducted,  and  the  results  compared  with 
those  achieved  using  holography. 

The  area  of  heat  transfer  is  relatively  untouched  by  the 
literature  in  conjunction  with  the  rise  of  speckle  metrology. 
It  is  felt  that  transient,  two-dimensional,  heat  transfer 
problems  could  be  analyzed  using  speckle  metrology  methods. 

Analyses  of  transient  and  real  time  phenomena  using 
television  speckle  recordings  have  been  conducted  with 


relative  success  by  both  Maddux  /” 1_7  and  German  Da  Costa 
of  Simon  Bolivar  University  /~21_7»  This  area  of  study  could 
be  undertaken  for  less  than  $5,000,  greatly  expanding  the 
field  of  speckle  metrology  for  the  Naval  Postgraduate  School. 

Currently  in  industry  speckle  metrology  is  being  utilized 
for  crack  detection  and  non-destructive  testing  /” 21_7, 
pressure  vessel  inspection  /~ 2 1_7,  and  stress  concentration 
analysis  /”21_7.  The  applicability  of  speckle  metrology  to 
Naval  Engineering  should  be  highly  evident  by  the  remarks 
entertained  in  these  recommendations  and  conclusions.  Again, 
as  mentioned  previously,  this  is  a  relatively  new  science; 
the  potential  for  growth  is  constrained  by  the  intensity  of 
the  analysts'  inquisitiveness.  Further  studies  are  con¬ 
sidered  a  must,  and  they  are  assumed  to  be  very  enlightening 
and  positive  contributors  to  modern  science  and  Naval  Engineering. 


APPENDIX  A 


Notes  Added  During  Proof 

After  the  thesis  was  approved ,  the  author  continued 
research  because  of  extreme  interest  in  the  development  of 
speckle  metrology.  The  image  shearing  camera  described  in 
the  Strain  Analysis  section  of  this  thesis  was  used  to  ob¬ 
tain  slopes  of  modal  amplitudes  at  various  frequencies  and 
compare  these  results  with  those  achieved  by  Fahey  /A-l7. 

I .  VIBRATION  THEORY 

The  specimen  of  Fig.  41  undergoes  a  steady  state 
sinusoidal  oscillation  represented  by: 

w(x,y,t)  »  A(x,y)  cos  (cut  +  0) 
where  w(x,yft)  ■  normal  displacement  of  a  point  P(xry)  at 

a  time  t. 

A(x,y)  »  modal  amplitude  function 

<u  ■  circular  frequency 

6  ■  arbitrary  phase 

and  if  normal  incidence  of  light  is  employed,  then  A 
is  a  function  of  time  given  by: 

l(x,y,tl  -  - 

Thus  IjfXty)  is  a  time  dependent  function  given  by, 

-  2A2  [l  +  cos2  |28x'  +  f^xSy') 

-  f2  (x’,y')  +  A(x',y',t)]] 
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where  I2  is  the  intensity  distribution  received  by  the 
photographic  plate.  If  the  film  is  exposed  to  I^Cx'^y'^t) 
for  a  time  T  which  is  large  compared  to  the  period  of 
vibration/  the  recorded  intensity  is  an  integration  of 
I2(x,/y'/t)  over  the  time  T  given  by: 

lifx'/y')  -  /0T  I2(x\y*,t)dt 


where  I^Cx'/Y1)  is  the  integrated  intensity.  The  above 
integration  yields: 


Ii(x'/y')  -  T  [l  +  cosy—  2Bx*  +  f^x'/y')  -  f2(x'/y*)J 

r  3ir  »„  3A(x#y)  1 

J0  Ll“  5x— 3x—  J 

This  is  a  case  where  a  high  frequency  term  is  amplitude 
modulated  by  a  zero  order  Bessel  function. 

Thus  Moire  fringes  will  be  formed  when: 

Jo  [f*«»  u,x'y)  ]  *° 


Again/  the  carrier  fringes  are  not  visible  by  the  eye. 
Therefore  Fourier  filtering  is  required  to  make  the  fringes 
visible.  See  Fig.  39. 

Figure  41  shows  the  laboratory  equipment  arrangement 
for  imaging  the  slopes  of  the  modal  amplitudes.  The  equip¬ 
ment  arrangement  is  identical  to  the  setup  of  Fig.  38  for 
strain  analysis  with  one  exception.  A  piezoelectric  excitor 
was  placed  behind  the  plate  and  excited  the  center  of  the 
plate  at  various  frequencies  and  amplitudes  using  the  oscil¬ 
lator  and  transformer  described  previously  in  the  vibrational 
section  of  this  thesis. 


140 


Figure  41.  Laboratory  equipment  arrangement  for  vibrational  analysis  using  the 
image  shearing  camera. 


Calculations  were  made  previously  by  Fahey  /”A-l7 
from  Leissa  £”a-2 7  for  the  frequency  parameters  fR  for  an 
all  sides  clamped  rectangular  plate. 
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The  fundamental  frequencies  calculated  by  Fahey  /~ A- 17 
were: 


“1,1 

9 

327.0  Hz 

“1,2 

9 

506.2  Hz 

“2,1 

- 

794.8  Hz 

“2,2 

« 

968.45  Hz 

Fundamental  frequencies  were  found  experimentally  to  be: 


“i.i 

m 

258 

Hz 

“1.2 
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548 

HZ 

H 

CM 

3 

- 

752 

HZ 

“2,2 

- 

1047 

HZ 

Modal  amplitude  fringes  were  obtained  for  the  three  highest 
frequency  responses,  while  at  258  Hz  fringes  were  unable  to 
be  obtained  due  to  the  low  amplitude  of  vibration  achieved 
by  the  excitor.  (See  Fig.  42a).  A  dark  ring  is  observed  at 
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(c)  Slope  of  modal  amplitude  at  752  Hz. 


(d)  Slope  of  modal  amplitude  at  1047  Hz. 


Figure  42  (Continued) 
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outer  plate  where  the  circular  first  mode  fringes  would  have 
been  seen  had  a  larger  amplitude  been  achieved. 

As  observed  previously/  although  the  fringes  produced  by 
holography  are  much  clearer  no  vibrational  isolation  is 
needed  nor  is  the  precision  optical  setup  required  for  speckle 
photography.  A  relatively  large  range  of  amplitude  variation 
was  achieved  using  the  image  shearing  camera  technique. 

Further  research  is  definitely  required  to  enable  the 
author  to  make  additional  conclusions  as  to  the  effectiveness 
of  this  technique. 
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